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PREFACE 


This  paper  constitutes  a departure  from  the  previous,  formal  and 
analytical,  treatments  of  the  subject.  It  presents  little  that  is  novel; 
it  is  written  colloquially;  it  is  concerned  just  as  much  with  mere  verbal 
qualitative  discussions  as  with  quantitative  evaluations;  it  is  frankly 
concerned  more  with  reader's  attitude  to  the  facts  than  with  the  facts 
themselves;  and  it  makes  no  pretense  of  brevity.  Nevertheless,  this 
presentation  has  been  found  necessary.  There  are  many  concepts  which  ought 
to  be  more  widely  understood  if  our  thinking  is  not  to  be  ossified,  and  if  we 
are  to  continue  - ty  working  together  - to  push  out  the  frontiers  of  ballistics. 
In  the  present  state  of  art  not  all  of  these  concepts  can  be  put  in  the 
precise  language  of  mathematics;  but  to  continue  to  ignore  them  would  be  a 
poor  substitute  for  precision.  Verbal  discussion  and  studied  colloquialism 
have  a legitimate  - if  modest  - place  in  ballistics.  Certainly,  facts  are  of 
no  use  unless  they  are  known  to  those  who  have  an  opportunity  to  use  them. 
Brevity  is  often  overdone:  a ballistician  spends  many  of  his  evenings  with 
pencil  and  paper,  trying  to  understand  the  unstated  implications  of  a brief  — 
and  therefore  allegedly  easy  - text. 

The  premise  of  this  paper  is  that  a visualisation  of  a physical 
phenomenon  is  useful:  not  so  much  because  on  occasions  it  might  be  all  that 
is  necessary,  as  because  on  few  occasions  it  might  open  a new  vista.  Of 
course,  visualization  is  part -arid-parcel  of  ary  understanding;  but  the  reader 
kncws  himself  that  it  is  all  too  often  crowded  out  by  algebra. 

The  paper  is  meant  to  supplement,  rather  than  to  supplant,  the  existing 
texts  on  the  subject.  Thus,  the  avoidance  of  the  calculus  is,  of  course, 
merely  a pose;  but  it  is  a pose  which  does  noi  harm,  and  which  might-  on 
occasions  - be  stimulating. 
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A SIMPLIFIED  APffiOACH  TO  THE  WING  MOTION  OF  A SPINNING  SHELL 


ABSTRACT 

An  attempt  is  made  to  free  the  theory  of  yawing  motion  of  a shell 
from  the  excessive  and  frustrating  reliance  upon  the  assumptions  of 
linearity.  A method  of  looking  upon  the  yawing  motion  is  suggested 
whereby  the  physical  significance,  or  the  mechanism,  of  the  effect 
produced  by  various  forces  and  torques  acting  upon  the  shell  can  be 
vividly  visualized,  without  resorting  to  the  linearization  of  the 
problem  (or  to  calculus). 

The  method  consists  essentially  of  considering  the  vector  of 
angnlar  momentum  L,  and  separating  the  dynamic  concept  (the  motion  of 
L under  the  influence  of  applied  torque,  here  termed  “quasi-precession® ) 
From  the  kinematic  concept  (the  motion  of  the  shell  with  respect  to  L, . 
here  termed  "quasi-outation®). 

The  method  is  tested  by  an  application  to  the  solved,  linear,  case: 
it  is  used  to  derive  - and  to  interpret  - all  conclusions  of  the  linear 
theory. 

The  principal  value  of  the  method  is  pedagogical.  It  shows  a promise 
in  tackling  certain  more  complicated  problems  of  ballistics,  such  as 
stability  of  liquid-filled  shell. 
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Unit  vector  along  shell  axis;  A « axial  moment  of  inertia 

Moment  of  inertia  about  a transverse  axis  through  the  c.g*  and 
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Vector  of  the  sum  of  all  applied  torques* 
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Vector  of  angular  momentum* 

Mass  of  body* 
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The  normal  force* 

The  c*g* 

Quasi-nutation  (the  angle  A^oL) * 

Quasi-precession  (the  angle  LoT^)* 

Synge's  force* 

Linear  Stability  factor* 

Unit  vector  along  trajectory;.  T ■ nagnus  torque 
Velocity  of  c*g* 

An  auxiliary  angle  as  shown  on  sketches* 

An  auxiliary  angle  as  shown  on  sketches* 

The  yaw. 

The  quasi-precession* 

Component  of  spin  (a)  perpendicular  to  longitudinal  axis  of  body* 
G>ef£ icient  associated  with  M,  define  $i  by  M ■ p-sinfi  or  M ■ /a  6, 

Taw  (viewed  as  a complex  number  in  the  linear  theory) * 

The  rate  of  change  of  yaw* 

Density  of  the  air* 
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Quasi-nutation* 
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INTRODUCTION 


The  theory  of  the  yawing  motion  of  a shell  is  an  ancient  and  important 
subject  of  ballistics s it  lies  at  the  root  of  well-nigh  all  problems  of 
ordnance  that  involve  the  accuracy  of  fire.  It  seems  that  this  theory 
deserves  being  known  more  widely,  and  used  more  frequently,  than  is  the 
case  now.  Also,  it  deserves  an  occasional  re-assessment. 

Traditionally,  the  subject  has  been  the  domain  - and  perhaps,  monopoly  - 
of  what  we  shall  call  the  mathematical  approach.  This  situation,  apparently, 
is  the  result  of  the  elegance  and  power  with  which  this  approach  solves  the 
basic  problems  of  the  stability  of  a shell.  ¥e  now  take  it  far  granted  that 
the  first  step  in  any  investigation  of  the  yawing  motion  must  be  an  idealization 
of  the  problems  that  the  elementary  physical  facts  must  be  condensed,  in  a 
simplified  form,  into  the  mathematical  notation;  that  the  interplays  of  these 
facts  can  be  handled  as  mere  algebraic  substitutions,  and  the  totality  of 
these  facts  must  be  expressed  as  a system  of  differential  equations,  the 
solution  of  which  may  be  obtained  by  abstract  mathematical  methods,  nob 
requiring  an  initmate  understanding  of  the  physical  significance  of  each 
intermediate  step;  that  the  physical  significance,  if  and  when  needed, will 
re-emerge  in  examining  the  mathematical  solution;  that  the  experimental 
tests  must  be  designed,  and  the  experimental  data  examined,  around  this 
abstract  solution;  and  finally,  that  predictions  may  be  made  by  these  methods 
that  will  possess  considerable  generality. 

It  is  of  interest  to  observe  that  this  tradition  appears  to  be  relatively 
young;  for  lb  is  difficult  to  trace  it  far  beyond  the  classical  work  of  Fowler 
in  the  First  World  War.  Raised  in  the  belief  in  the  power  of  mathematical 
deduction,  we  might  surmise,  far  instance,  that  the  introduction  of  rifling 
was  a result  of  Ruler's  theory  of  spinning  top;  this  indeed  might  have  been 
relevant  as  far  as  artillery  is  concerned.  But  the  rifles  of  the  American 
Revolution  were  made  by  the  village  blacksmith,  innocent  of  the  theory,  while 
the  most  advanced  nation  in  Europe  was  still  mass-producing  the  smooth-bore 
Brown  Bess.  The  importance  of  the  mathematical  approach,  clearly,  rests  not 
upon  tradition,  but  on  hopes  that  it  will  produce  useful  results  sooner  than 
a haphazard  invention  and  empiricism  can  produce  them.  It  is  therefore  not 
inappropriate  that  from  the  proving-ground  point  of  view  ("the  proof  of  the 
pudding  is  in  the  eating")  we  inquire,  occasionally,  whether  the  mathematical 
approach  itself  - that  is,  the  idea  that  it  is  possible  to  apply  the  rigor 
of  mathematics  to  a field  of  engineering  (to  wit,  ballistics;  - has  lived 
up  to  the  expectations. 

The  fact  is,  the  problem  of  the  yawing  motion  of  a shell  has  been  solved 
(and  indeed,  vindicated),  only  for  the  highly  simplified,  linear  case.  Since 
the  underlying  assumption  usually  is  that  the  yaw  is  small,  this  is  as  though 
the  theory  solves  the  problem  only  when  the  shell  flies  rather  satisfactorily. 

When  some  of  the  simplifying  assumptions  (such  as  the  linearity  of  the 
force  and  torque  system,  smallness  of  the  yaw,  symmetry  of  the  shell,  etc. ) 
have  to  be  rescinded,  the  complexity  of  the  problem  grows  rapidly.  In  fact, 
most  of  the  present  work  in  this  field  anticipates  an  extensive  use  of 
modern  high-speed  computing  machines.  In  this  connection  it  is  particularly 
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annoying  that  with  the  recis&on  of  the  linearity  it  in  no  longer 
legitimate  to  utilize  that  basic  feature  of  the  linear  theory,  the  de- 
composition of  the  yaw  into  two  normal  modes,  precession  and  natation* 

The  numerical  solution  of  the  equations  of  yawing  motion  becomes  extremely 
laborious,  and  yet  much  of  the  information  expected  from  this  solution 
(e.g.,  the  information  that  the  shell  is  yawing  approximately  epicyjLoidaliy) 
resembles  very  much  the  information  already  available  from  the  linear  theory. 
The  practical  result  usually  is  that  the  available  indications  of  nan-linearity 
are  simply  suppressed  - rather  than  investigated.  In  other  words  (to  mix 
a metaphor),  the  elegance  of  the  linear  theory  has  exerted  such  influence 
on  the  methods  of  ballistics  that  the  present-day  theory  is  chained  to  the 
assumptions  of  linearity  as  to  a point  of  diminishing  returns.  There  exists 
an  acute  need  far  a more  powerful  method  of  tackling  the  problem  of  yawing 
motion. 

The  method  suggested  hereby  has  no  claim  of  being  a complete  answer 
to  this  need;  yet  it  might  he  a contributing  factor.  What  it  has  to  offer 
is  simply  an  easier  visualization  of  the  problem ; and  it  is  quite  possible 
that  a future  method  will  be  dependent  upon  a visualization  and  a judgment*. 

In  one  problem,  particularly  - that  of  liquid-filled  shell  - the 
traditional  mathematical  approach  appears  to  be  simply  stalemated.  It  has 
not  been  possible  to  get  even  to  first  base,  i.e.,  to  formulate  total 
differential  equations  of  the  motion;  while  a solution  of  partial  differential 
equations  of  the  motion  of  the  liquid  within  the  shell  (with  the  object  of 
using  the  results  in  determining  the  motion  of  the  shell)  seems  a long  way 
off.  It  is  in  connection  with  this  problem  that  the  method  suggested  in  this 
report  arose.  It  started  as  an  attempt  to  by-pass  the  differential  equations 
by  means  of  visualizing  the  motion  of  the  liquid  within  the  yawing  shell, 
or  as  an  esthetic,  or  qualitative  "reasoning11  (if  such  a term  may  be  applied). 
It  was  not  easy,  however,  to  meet  the  contemporary  standards  set  by  the 
tradition  of  the  mathematical  approach;  in  fact,  a number  of  pitfalls  of 
this  "qualitative  reasoning”  have  been  found.  Far  this  reason  it  appeared 
best  to  expound  this  method  far  the  simple  case  of  a rigid  shell,  and  to 
test  it  by  applying  it  to  the  solved  linear  case  (where  this  method  can 
be  used  to  re -derive  and  to  interpret  all  conclusions  of  the  linear  theory); 
postponing  the  possible  application  of  this  method  to  the  cases  of 
non— linearity,  asymmetry,  and  liquid  filler. 


FHELIMIMRIES 

Our  method  is  non -mathematical  only  in  this  respect;  we  shall  not 
require  of  the  reader  the  formal  knowledge  of  calculus  and  of  differential 
equations.  This;  Is  rather  a tenuous  simplification,  since  the  differential 
equations,  of  course,  are  the  essence  of  the  problem.  On  the  other  hand, 
we  shall  expect  the  reader  to  have  a mastery  of  algebra,  geometry,  plane 
(and  possibly,  spherical)  trigonometry;  of  the  concepts  of  torque  and 
angular  velocity  as  vectors  defined  in  the  right-handed  sense;  of  the 
principal  axes  of  inertia  (i.e.,  the  understanding  that  a rigid  body  may 
spin  steadily  only  if  it  is  dynamically  balanced  about  the  axis,  and  that 
there  are  at  least  three  such,  mutually  perpendicular,  axes  in  every  rigid 
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body)  | of  the  moments  of  inertia;  of  the  concept  of  the  vector  of  angular 
momentum  (which  can  be  constructed  by  decomposing  the  instantaneous  vector 
of  angular  velocity  along  the  three  mutually  perpendicular  axes  of  inertia 
of  the  body,  multiplying  each  component  by  the  corresponding  moment  of 
inertia,  and  recombining  the  products  vectorially)  j of  the  fact  that  the 
vector  of  applied  torque  is  the  velocity  of  the  tip  of  the  vector  of  angular 
momentum;  of  the  resolution  of  the  motion  of  a rigid  body  into  the  motion 
of  its  center  of  gravity  and  the  rotation  about  this  center  of  gravity; 
and  of  the  use  of  complex  numbers. 

We  shall  consider  only  an  axially -symmetr ic , spindle-like  body  with 
moments  of  inertia  A and  B (A  B,  B the  same  about  every  transverse  axis 
through  the  center  of  gravity).  For  such  a body  the  unit  vector  along  the 
axis  (A,),  the  vector  of  angular  velocity  (©)  and  the  vector  of  angular 
momentum  (1)  are  always  co-planar  (Figure  ij*  For,  if  the  instantaneous 
motion  of  the  body  (with  respect  to  its  center  cf  gravity)  is  a rotation 
at  an  angular  velocity  © about  an  axis  (to)  that  is  different  from  A^,  the 

vector  © can  be  resolved  into  the  axial  spin  along  A^  and  the  cross-spin 

Tj  perpendicular  to  A^;  because  cf  the  symmetry,  the  direction  ^ of  vj  may  be 

taken  as  the  second  principal  axis  of  inertia  (there  being  no  component  of  © 
along  the  third  axis),  and  L can  be  constructed  from  the  component  Aco^  along 

A^  and  along  7^  , i,e,,  in  the  plane  of  and  co,  In  particular,  since 

B^A,  L is  always  further  away  from  than  a is,  and  there  exists  the  relation 


tan(A1,  a)  rj  /c^ 

l)  * ' k/B 


Note  that  ,^his  relation  allows  a coincidence  of  the  vectors  © and  L 
(which  is  to  say,  allows  a steady  unrestrained  rotation  of  the  bo3y,  or“a 
steady  rotation  with  no  torque  applied)  only  when  the  angle  between  A,  and 

» ie  0,  n or  n/2*  in  which  cases  the  motion  is  either  a pure  spin  (right- 
handed  or  left-handed)  or  a pure  cartwheeling. 

Generally,  however,  the  vector  of  angular  velocity  <0  is  fixed  neither 
in  space  nor  in  the  body.  It  is  for  this  reason,  perhaps,  that  students 
of  ballistics  have  same  slight  difficulties  in  visualizing  the  three- 
dimensional  yawing  motion  of  a rigid  body. 

The  visualization  is  facilitated  by  the  decomposition  of  the  vector 
© into  the  components  and  TJ  , The  cross-spin  7)  , of  course,  is  simply 

the  velocity  of  the  "tip  of  the  unit  vector  A^  in  the  surface  of  the  sphere 

which  is  centered  at  the  c,g,  of  the  body,  and  moves  with  the  body,  but 
does  not  rotate;  this  velocity,  of  course,  is  at  right  angles  to  the  vector 
7)  • In  this  way  we  replace  the  visualization  of  the  motion  of  the  body 
{rotation  about  ©)  by  a visualization  of  the  motion  of  a line  A. , with  the 
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Fig.  I:  Angular  Velocity  4 Angular  Momentum 
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axial  spin  superposed  if  and  when  necessary.  In  other  words,  a three- 

degree-of -freedom  situation  is  replaced  by  a two-dimensional  situation 
(in  the  surface  of  the  unit  sphere)#  The  motion  of  the  tip  of  the  unit 
vector  in  the  surface  of  the  unit  sphere  can  be  specified,  of  course, 

as  a rotation  about  the  vector  co. 


However,  the  vector  <a  is  not  a convenient  point,  of  reference*  The 
vector  L is  much  more  convenient.  We  shall  presently  -see  that  it  does 
not  move  as  fast  as  co  and  in  particular,  we  knew  how  it  does  move  - viz., 
the  velocity  of  its  £ip  is  the  vector  G of  the  applied  torque*  The  motion 
of  the  tip  of  may  be  described  as  a- rotation  about  the  instantaneous 

position  of  1 at  an  angular  velocity 

•d  ■ a g * * 1/B*  (2) 

where  L is  the  magnitude  of  L* 

We  can  thus  describe  the  motion  of  the  body  ty  two  statements: 

I,  A^  rotates  about  the  instantaneous  position  of  L at  an 

angular  velocity^,  and 

U.  The  velocity  of  the  tip  of  L is  G, 

it  being  understood  that  G can  be  determined  from  the  specification  of 
A^  and  L (and  other  relevant  and  available  data).  The  statement  U,  of 

course,  is  simply  the  applicable  law  of  physics;  statement  X connects  the 
easily  visualizable  (and  measurable)  vector  A^  with  the  very  relevant 

physical  quantity,  L.  In  most  texts  on  ballistics  the  vector  of  angular 
momentum  is  mentioned  briefly  in  the  beginning  of  a derivation,  and 
thereafter  is  replaced  by  the  appropriate  combinations  of  its  elements 
(A,  B,  and  T)  ).  Our  suggestion  amounts  to  inviting  the  reader  to 

* think  in  terms  of  angular  momentum",  i*e.,  to  make  a better  - mar© 
frequent,  end  more  intimate  - use  of  this  essential  oonoept.  The  two 
statements,  together  with  a specification  of  G correspond,  of  course, 
to  the  differential  equations  of  the  angular  motion  of  the  shell;  they 
are  to  be  augmented  by  the  statement  that  the  acceleration  of  the  center 
of  gravity  is  the  applied  force  divided  by  the  mass  (with  force  properly 
specified; • 


VACUUM  NUTATION 

An  hearten t special  case  occurs  when  there  are  no  torques  acting 
upon  the  body  (e.g*,  shell  flies  in  vacuum)*  The  the  vector  L is 
constant.  A very  simple  description  of  the  motion  is  then  possible. 

As  the  vector  on  Figure  1 comes  out  of  the  plane  of  the  paper  (rotating 

about  a>),  the  co-planarity  of  A. , to  and  L must  be  preserved;  hence  the 
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vector  os  coir.es  out  of  the  plane  of  the  paper*  too.  Since  the  angle  A-^OL 

is  not  changing,  by  (1)  the  angles  A^Qgj  and  aOL  are  not  changing  either; 

hence  the  whole  triplet  of  vectors  A^,  co  and  L ia  rotating  (as  a rigid 

plane)  about  L»  The  motion  of  the  body  is  a rolling  of  the  circular  cone 
of  haif-angle*~AjOca  on  the  circular  cone  of  half -angle  coOI*. 

It  is  not  amiss  at  this  point  to  make  several  trivial  observations. 
The  motion  is  not  a simple  rotation,  it  is  rolling;  yet  it  is  a basic, 
steady,  type  of  motion,  that  may  go  cut  forever  - if  our  assumptions  (so 
far,  only  the  rigidity  and  symmetry  of  the  body  and  the  absence  of  applied 
torque)  hold;  this  is  to  say  that  the  motion  can  be  predicted  for  any 
length  of  time.  If  the  assumptions  are  not  quite  exact,  the  motion  still 
can  be  predicted,  with  seme  accuracy,  far  some  interval  of  time  — the 
more  nearly  true  these  assumptions,  the  greater  this  interval.  The  angle 
AjOjLnay  have  any  value  from  0 to  n;  it  may  also  have  an  arbitrary  initial 
orientation  (i.e.,  there  are  two  scalars  representing  the  arbitrary 
constants  of  integration). 

We  shall  call  the  angle  A^OL  B quasi-nut at ion" . 

Since  in  this  problem  the  only  direction  fixed  in  space  is  L,  viz., 
we  may  visualize  the  shell's  c.  g.  as  stationary  in  space, it  might  be  also 
natural  to  call  this  angle  - in  this  case  - the  yaw.  Thus  the  quasi- 
nutation could  be  defined  as  the  yaw  with  respect  to  Ij  or,  we  might  say 
that  in  vacuum  nutation  the  yaw  with  respect  to  L consists  of  quasi- 
nutation. ““ 


OVERTURNING  MOMENT 

The  torque  (or  moment)  acting  upon  the  shell  in  flight  is  specified, 
conventionally,  as  the  sum  of  its  components;  these  components,  in  turn, 
being  specified  both  in  direction  and  magnitude,  in  certain  convenient  ways. 
It  is  not  quite  a simple  matter  to  outline,  at  the  outsat,  which  components 
(or,  as  the  parlance  goes,  which  torques)  should  be  considered;  the  student 
of  ballistics  would  find  it  best  to  go  by  easy  steps,  digesting  well  the 
basic  concepts,  and  accepting  or  rejecting  the  more  advanced  concepts  as 
the  need  for  them,  as  well  as  the  practicability  of  handling  them,  becomes 
clear. 

Of  these  torques  the  most  important  one,  theoretically,  is  the  one 
whose  vector  is  perpendicular  to  the  plane  of  yaw  (i.e.,  whose  couple 
is  in  the  plane  of  yaw),  and  whose  magnitude  depends  essentially  upon 
the  angle  of  yaw  (angle  from  trajectory  of  center  of  gravity  to  axis  of 
shell).  It  is  called  usually  simply  "moment*1,  since  historically  it 
was  the  only  moment  at  first  considered;  and,  since  it  was  first  considered 
in  connection  with  spin-stabilization  of  artillery  shell,  it  is  reckoned, 
conventionally,  positive  when  overturning.  Hence  the  mare  exact  name 
(distinguishing  it  from  other  components;  is  overturning  torque.  With 
fin-stabilized  projectiles  the  righting  torque  is  reckoned  merely  as  a 
negative  overturning  torque. 
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In  the  linear  theory  the  magnitude  of  this  torque  (M)  is  assumed  to 
be  proportional  to  the  angle  of  yaw  (6).  The  solution  then  is  remarkably 
elegant*  In  certain  texts  on  ballistics  M is  assumed  to  be  (for  the 
reasons  of  certain  convenience  in  the  preliminary  mathematics)  proportional 
to  sin6  (often  a misleading  assumption)  or  tan6  (sometimes  a better 
approximation,  though  still  limited  to  small  angles  of  yaw) 5 in  practice, 
both  of  these  assumptions  eventually  revert  to  the  assumption  of  linearity 
in  6.  The  exact  (and  often  not  simple)  function  M(6)  can  be  obtained 
experimentally,  say,  ty  wind-tunnel  tests,  in  exact  mathematical  solution 
has  been  worked  out  but  is  - in  comparison  with  the  linear  theory  - so 
cumbersome  as  to  be  considered  impractical.  It  is  precisely  this  gap 
between  the  linear  and  non-linear  theories  which  our  approach  attempts 
to  bridge  (or  more  exactly,  our  approach  attempts  to  outline  a method 
of  bridging  this  gap). 


MUTATION  AND  PRECESSION 

There  are  two  special  cases  when  the  yawing  motion  Is  particularly 
simple  - to  wit,  circular  - even  with  a non-linear  M(6).  These  cases 
arise  particularly  naturally  if  we  assume,  for  the  time  being,  that  K 
is  very  small,  i.e.,  1 moves  only  slowly;  the  motion  is  almost,  but  not 
quite,  the  same  as  the  nutation  in  vacuum. 

Let  T ^ be  a unit  vector  in  the  direction  of  the  trajectory,  i.e., 
in  the  direction  of  the  velocity  u of  the  center  of  gravity  of  the  shell. 

Consider  first  the  case  when  T^  is  approximately  in  the  same  direction 

as  X.  The  situation  may  be  such  (Figures  2a  and  2b)  that  the  triplet  of 
vectors  A^,  © and  1 rotates,  as  a rigid  plane,  about  T^.  In  particular, 

if  M is  an  overturning  torque  (Figure  2a),  L would  be  coming  out  of  the 
plane  of  the  paper;  the  motion  then  may  be  steady  if  T^  is  farther  from 

than  1 is.  Or,  if  M is  a righting  torque  (Figure  2b),  X would  be 

moving  into  the  plane  of  the  paper;  but  the  motion  still  may  be  steady 
if  is  closer  to  A^  than  X is.  In  either  case  the  notion  of  the  body 

is  a rolling  of  thd  cone  A^Cto  (same  as  before),  but  this  time  not  on  the 

cone  of  half -angle  ©01  (as  was  the  case  in  vacuum)  - rather,  on  the  cone 
of  half -angle  coOI^.  The  angular  velocity  of  the  plane  of  yaw  is  slightly 

lower  than  *fl  with  spin-stabilized  shell  (Figure  2a),  slightly  higher  than 
.flwith  fin-stabilized  shell  (Figure  2b). 

Such  a motion  is  call  nutation.  In  it  the  yaw  (the  angle  A^OT^)  is 

approximately  the  same  as  quasi-nutation  (angle  A^QL),  but  also  has  a 

slight  component  T^QL.  We  shall  call  the  angle  T^OL  "quasi-precession11. 

Note  that  with  both  spin-stabilized  and  fin-stabilized  shell  the  nutation 
is  in  the  same  direction  as  the  spin. 


15 


Notation  • O denotes  the  tip  of  arrow  coming  out  of  the  paper 

® denotes  the  tail  (feathers)  of  arrow  going  into  the  paper 


spin  - stabilized 


shell  fin -stabilized  shell 

(Ti  between  A,  $ L) 

Fig.  2:  Nutation 


spi n - stabi  H zed  shell 

Fig.  3 : 


fin  -stabi  l i ze  d shell 

(T,  outside  of  angle  AQL ; left-handed  yawing) 

Precession 


Consider  next  the  case  when  is  far  away  from  the  triplet  of  vectors 

A.,  co  and  L.  This  might  be  the  case,  for  instance,  when  the  angle  A^OL  is 

small'  and  this  would  mean  that  the  cross-spin  T)  is  small,  i.e.,  the  shell  axis 
A^  is  moving  only  slowly.  The  situation  again  may  be  such  (Figures  3a  and 

3b)  that  the  triplet  A^,  co  and  L rotates  as  a rigid  plane  about  T^.  In 

particular,  if  M is  an  overturning  torque,  T^  is  still  on  the  same  side  of 

as  L is,  though  much  farther  away  (Figure  3a)#  Or,  if  H is  a righting 

torque,  T^  must  be  on  the  other  side  of  from  L (Figure  3b).  The  motion, 

again,  is  a rolling  of  the  same  cone  A^Oco  on  the  cone  coOT^.  For  spin- 

stabilized  shell  the  angular  velocity  of  the  plane  of  yaw  is  in  the  same 
direction  as  spin,  but  much  slower  than  the  rate  of  vacuum  nutation;  for 
fin-stabilized  shell,  however,  the  motion  is  in  a direction  opposite  to 
spin. 


Such  a motion  is  called  precession.  In  it  the  yaw  consists  mostly 
of  the  quasi-precession,  but  there  is  also  a slight  component  of  quasi- 
nutation. 


Nutation  and  precession  are  obviously  pretty  much  the  same  kind  of 
phenomenon;  the  difference  between  them  appears  to  be  merely  quantitative. 
In  fact,  this  seems  to  be  precisely  the  case  as  far  as  the  original, 
astronomical,  usage  of  these  wards  is  concerned;  the  phenomena  which 
go  by  these  names  are  of  the  same  nature,  viz.,  the  substantially 
circular  yawing  motion  of  an  axis  of  a spinning  symmetrical  rigid  body 
under  the  influence  of  certain  external  torques.  Yet  it  would  not  be 
right  to  ignore  the  sense  of  the  qualitative  distinction,  conveyed 
by  the  very  choice  of  the  words?  precession,  coming  from  the  sane  root 
as  "precede11,  "proceed",  "procession",  emphasizes  the  stateliness,  viz., 
steadiness  and  slowness,  of  the  motion;  while  natation  (meaning  nodding) 
emphasizes  the  disturbance,  the  unsteadiness,  - i.e#.  perturbation,  or 
oscillation.  Thus,  while  we  now  find  it  convenient  (from  the  mathematical 
point  of  view)  to  consider  these  two  phenomena  separately,  and  as  being 
on  a par,  it  is  essential  to  keep  in  mind  that  historically  they  arose 
as  components  of  a larger,  more  complicated,  motion;  and  that  the  pre- 
cession is  somehow  more  fundamental,  the  nutation  being  of  the  nature  of 
a refinement  on  the  precession. 

In  our  case  the  qualitative  distinction  is  that  the  nutation 
emphasizes  the  motion  of  shell  axis  about  the  vector  of  angular  momentum 
(so  to  say,  a kinematic  concept),  while  precession  emphasizes  the  motion 
of  the  vector  of  angular  momentum  about  some  other  axis,  in  this  case, 
trajectory  (a  dynamic  concept).  Our  concepts  of  "quasi-nutation"  and 
"quasi-precegsion"  serve  further  to  crystallize  this  qualitative 
distinction. 


These  terms  should  be  acceptable  to  purists?  "quasi®  means  "something 
like",  and  besides,  like  "nutation"  and  precession", -is  Latin;  while 
"pseudo"  means  "falsely  appearing  as",  and  is  Greek.  Yet  the  terms  are 
still  rather  vague?  quasi-nutation  is  really  meant  to  be  the  essence  of 
nutation,  and  quasi-precession,  the  essence  of  precession. 
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While  the  nutation  and  the  precession,  so  defined,  do  not  depend  on 
the  assumption  of  linearity  of  M(6),  neither  do  they  give  us  much  help  in 
passing  from  these  two  special  cases  to  the  general  (epicyclic)  case  of 
* 

yawing  motion.  Their  importance,  however,  is  very  much  enhanced  once 
the  assumption  of  linearity  of  M(b),  i.e.,  the  smallness  of  is  made. 

Then  their  angular  rates  (and  also,  as  the  mare  elaborate  theory  shows, 
their  rates  of  damping)  become  independent  of  the  angle  of  yew,  and  these 
two  types  of  motion  became  superposable.  Moreover,  they  then  become 
"normal  modes"  of  motion  (viz.,  not  only  independent  and  superposable, 
but  also  particularly  simple  modes).  In  fact,  the  elegance  of  the 
mathematical  approach  amounts  to  this:  in  the  linear  theory  it  is  possible 
to  decompose  any  complicated  (epicyclic)  motion  in  two  simple  components, 
nutation  and  precession. 

We  decompose  the  yaw  (the  angle  from  to  A^)  in  an  alternate  way, 

viz.,  into  quasi-precession  (the  angle  from  T,  to  L)  and  quasi-nutation 
(the  angle  from  L to  A^),  essentially  for  thereason  that  these  concepts 

survive  in  any  general,  non-linear,  theory*  i.e.  ^simply  bring  in  the 
vector  of  angular  momentum.  To  describe  the  yawing  motion,  we  must  follow 
the  motion  of  L (as  affected  by  the  yaw  and  other  factors),  superposing 
upon  it  the  rotation  of  the  quasi-nut at Ida? (about  the  instantaneous, 
moving,  position  of  L)  at  the  steady  angular  velocity  O* . (This  is 
equivalent  to  a numerical  solution  of  an  alternate  form  of  the  differential 
equations  of  motion).  This  decomposition,  however,  should  not  be  confused 
with  the  customary  decomposition  of  yaw  into  nutation  and  precession  (which 
is  possible  only  in  the  linear  theory).  The  quasi-nutation  and  quasi-preOession 
are  not  normal  modes  of  motion,  are  not  independent  of  each  other,  and  are  not 
constant  in  magnitude  and  rate  (as  the  nutation  and  precession  are).  In  the 
epicyclic  motion  the  quasi-nutation  (if  considered  separately)  and  the  quasi- 
precession  describe  not  the  simple  spirals,  bijfc  epicycles  (which,  however,  are 
closer  to  circles  that  the  epicycle  of  is  ).  An  additional  shortcoming 

of  our  presentation  is  that  for  slowly-spinning  f in-stabilized  projectiles 
the  magnitude  of  both  the  quasi-nutation  and  of  the  quasi-pre  cess  ion  is 
near  90°,  and  the  situation  is  somewhat  less  vivid,  since  the  (small)  yaw  ? 
is  effectively  a small  difference  between  two  large  quantities. 

While  we  spoke  of  these  quantities  (nutation,  precession,  quasi- 
nut  at  ion  and  quasi-precession),  rather  loosely,  as  "angles*,  thsy  are  not, 
of  course,  simple  angles  ( scalars) j ra^hery  each  one  of  them  should 


Unless  we  go  into  quite  advanced  mathematical  theory.  Recently  L.  H. 
Thomas  (Ballistic  Research  Laboratory  Report  No.  839,  Tfhe  Theory  of 
Spinning  Shell",  1 992)  has  introduced  - by  the  method  of  variation  of 
constants  - the  concepts  of  nutation  and  precession  into  the  general, 
non-linear,  theory. 
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properly  be  specified  by  two  numbers  (magnitude  and  orientation).  The 
situation  is  again  particularly  simple  in  the  linear  theory,  where 
these  quantities  are  complex  numbers,  or  "planar  "vectors",  and  are 
commutative  (the  order  in  which  they  are  added  is  immaterial).  Generally, 
these  concepts  are  not  very  elementary  (they  are,  strictly,  quaternions 
or  dyadics),  but  for  our  present  purposes  they  are  rather  obvious  concepts 
of  spherical  trigonometry. 


METHODS  QF  ffiESENTATlON 


The  interpretation  of  the  yawing  motion  in  terms  of  the  rolling  of 
one  cone  upon  another,  cannot,  unfortunately,  be  extended  to  the  general, 
non-circular,  yawing  motion.  We  therefore  have  to  resort  to  the  above- 
mentioned  abstraction,  viz.,  to  consider  the  motion  of  tbs  shell  axis 
A^  instead  of  the  motion  of  the  rigid  body  as  a whole.  The  latter 

motion  can  always  be  visualized  subsequently,  either  by  superposing  the 
axial  spin  upon  the  motion  of  A^,  or  by  constructing  the  vector  co 

from  the  instantaneous  positions  of  and  1,  and  (l). 


Let  us  imagine  a sphere  of  unit  radius,  centered  on  the  center  of 
gravity  of  the  shellj  the  sphere  moves  with  the  shell,  but  does  not 
rotate  ("remains  parallel  to  itself").  Let  us  view  the  interior  surface 
of  this  sphere  from  its  center,  and  let  us  attempt  to  sketch  the  relation- 
ship of  the  vectors  based  on  that  center,  as  it  would  then  appear.  The 
directions  in  space  would  now  map  as  points,  and  the  angle$  as  certain 
"straight"  lines  (the  arcs  of  great  circles).  On  such  a Sketch  we  will 
be  losing  the  distinction  between  the  vectors  of  certain  definite  length, 
and  the  unit  vectors  (which  terminate  on  the  surface  of  the  sphere,  and 
which  we  distinguish  by  the  subscript  l);  thus,  the  vector  of  the  velocity 
u of  the  shell  is  coincident  with  the  vector  T^,  etc.  The  exact  manner 


in  which  the  spherical  surface  is  mapped  on  our  sketches  is  immaterial 
at  this  stage:  the  surface  of  our  sketches  can  be  identified  with  the 
spherical  surface  only  when  the  angles  involved  are  small  (we  are  tten 
back  in  the  linear  theory),  and  whenever  the  angles  are  not  small,  we 
must  revert  to  the  t hr ee-Kiimens ional  representation.  Some  of  our  angles 
are  inherently  large:  e.g. , the  vector  M is  at  90°  to  both  T^  and  A^, 

vector  ^ is  at  90°  to  A^,  etc.}  it  would  be  awkward  to  represent  such 


vectors  as  points  on  our  plane  sketches,  but  we  can  readily  visualize 
them  as  arrows  projected  on  that  part  of  the  sphere  at  which  we  are 
looking.  Since  such  a projection  can  be  made,  so  to  say,  in  any- 
desired  direction,  the  vectors  so  projected  need  not  be  confined  to 
ary  particular  line  on  our  sketches,  but  must  have  only  a particular 
direction.* 


In  exact  terminology,  the  polar  vectors  (such  as  velocity  and  force) 
become  axial  vectors  (such  as  angular  velocity  and  torque)  on  our 
sketches.  Obviously,  we  shall  base  torques  or.  L,  forces  on  T^,  etc. 
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Let  us  look  forward  as  the  shell  flies,  i.e«,  let  us  center  our 
field  of  vision  on  the  point  (vector)  T^j  the  vector  does  not  have 

to  be  horizontal  (though  it  is  easiest  to  imagine  it  so).  If  the 
trajectory  is  a straight  line  (usually  not  a bad  assumption),  the  point 
T1  is  fixed  in  the  surface  of  our  sphere,  and  is  a very  convenient 

point  of  reference.  However,  if  the  trajectory  curves  or  swerves,  the 
point  moves.  The  velocity  of  is  easy  enough  to  compute*  The  tip 

of  the  vector  of  (linear)  momentum  of  the  shell,  (where  m is  the 

mass,  and  u the  velocity  of  the  shell),  moves  with  the  velocity  equal 
to  the  applied  total  (gravitational  and  aerodynamic)  force.  Hence, 
that  component  of  the  applied  force  which  is  perpendicular  to  T^, 

divided  by  mu,  is  the  (vectorial)  velocity  of  on  the  surface  of  our 

sphere.  Ve  should  also  bear  in  mind  that  the  total  force  may  have  a 
component  along  (e.g.,  drag)  which  can  cause  additional  effects. 

Similarly,  the  velocity  of  the  point  L (the  point  of  intersection 
of  vector  L with  the  surface  of  the  sphere)  is  the  component  of  the 
total  applied  torque  normal  to  L,  divided  by  L (the  magnitude  of  L); 
again,  there  may  exist  other  effects  due  to  a component  of  applied 
torque  along  L.  It  is  interesting  to  note  that  the  vector  of  applied 
torque  must  be  decomposed,  therefore,  along  and  across  L - while  it 
is  more  natural  (and  in  fact,  is  customary)  to  decompose  it  along  and 
across  A^. 

While  the  vector  of  cross-spin  'Tj  can  be  projected  as  an  arrow 
tangent  to  the  surface  of  the  sphere  at  the  point  A^,  it  is  simpler  to 

use,  instead,  an  arrow  representing  the  velocity  of  the  point  A^  in  the 

surface  of  the  sphere.  In  magnitude  this  is  '*)  s and  in  direction  this 
is,  of  course,  90°  counterclockwise  (on  our  sketches)  from  the  projection 
of  the  vector  T)  j and  90°  clockwise  from  LA^. 

In  order  to  systematize  somewhat  the  forbidding  possible  multiplicity 
of  the  arrangements  of  the  points  A^,  and  L on  our  sketches,  let  us 

also  adopt  this  conventions  let  us,  as  it  were,  tip  out  head  sidewise  in 
the  plane  of  yaw.  To  wit,  we  shall  draw  the  point  always  directly 

above  the  point  T^J  "this  will  not  mean  that  the  plane  of  yaw  is  vertical  - 

it  is  just  easier  to  imagine  it  so.  A naive  attempt  to  depict  the  reader 
in  the  process'  of  visualizing  an  instantaneous  motion  of  a shell  is 
given  in  Figure  U. 

Let  us  adopt  one  more  convention.  With  the  single  point  on  our 

sketches  representing  all  possible  yaws,  we  must  still  distinguish  various 
possible  directions  of  the  velocity  of  A,:  i.e.,  we  must  distinguish  the 

increasing  and  the  decreasing  yaws  (or,  an  instantaneously  circular  yawing 
motion);  as  well  as  right-handed,  left-handed  and  planar  yawing  motions. 
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definition 
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On  orur  sketches,  to  specify  the  direction  qualitatively  it  is  easiest 
to  refer  to  the  points  ctf  the  compass.  Let  us  then  repeat  each  sketch 
four  (or  eight)  times,  with  different  possible  directions  (N,  Mrf,  W,  Stf, 
S,  etc.)  of  the  velocity  of  A_;  and  let  us  place  each  component  sketch 
in  the  corresponding  direction  from  the  center  of  the  sketch  as  a whole . 
Accordingly,  our  sketches  will  be  arranged  in  the  following  scheme, 
whereby  the  instantaneous  yawing  motion  is  respectively: 


increasing 

left-handed 


Increasing 

instantaneously 

planar 


increasing 

right-handed 


instantaneously  circular 
left-handed 


instantaneously  circular 
right-handed 


decreasing 

left-handed 


decreasing 

instantaneously 

planar 


decreasing 

right-handed 


In  order  to  become  thoroughly  familiar  with  such  a method  of  presentation, 
let  us  new  consider  an  example:  let  us  investigate  qualitatively  the  mechanism 
of  the  familiar  looped  epicycle  of  a spin-stabilized  shell  with  right-handed 
spin* 


I®  CHAN  ISM  CF  AN  EPICXCLE 

Let  us  refer  to  Figure  5a,  and  in  particular,  commence  at  the  instant 
when  the  instantaneous  yawing  motion  is  circular  and  left-handed*  This 
corresponds  to  the  W sketch,  and  - as  it  will  be  presently  seen  on  Figure 
5 b - to-  the  minimum  of  yaw* 

Point  L now  must  be  above  A^*  The  overturning  torque  M moves  L to 

the  right,  and  the  situation,  obviously,  passes  into  that  described  on 
the  Ntf  sketch* 

Now  the  situation  has  very  pronounced  features  of  what  we  shall 
call  K local  instability”.  The  yaw  (T^A^)  is  growing,  and  particularly, 

in  growing  at  an  increasing  rate:  for,  not  only  is  L to  the  right  of 
A^,  but  also:  the  velocity  of  is  becoming,  on  cur  sketch,  mare 

nearly  vertical;  both  the  quasi-precession  and  - particularly  - quasi- 
nutation  are  growing,  as  though  stretched,  to  some  extent,  by  the 
torque  M,  so  that  with  the  constant  angular  velocity of  the  quasi - 
nutation  the  linear  velocity  of  is  Increasing;  and  finally,  as  the 

yaw  is  growing,  the  overturning  torque  M is  growing,  and  also  is  turning 
to  the  left,  and  is  thereby  speeding  up  the  velocity  of  L (it  might  be 
also  noted  that  presently  M begins  to  stretch  L,  increasing  thereby  the 
angular  velocity  <Clc£ o equation  (2)=  If  such  a situation  were  to  continue 
long,  the  yaw  would  be  quickly  built  up  to  very  high  values  (which  it 
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(h)  CL7?  example  : usual  epicycle 

Fig  5 : Method  of  Presentation 


Indeed  would  if  the  spin  were  too  lew.  So  that  L/B  is  small  and 
the  velocity  M/L  of  L is  large).  In  practice,  however,  the  situation 
(because  of  sufficiently  high  spin)  rapidly  passes  into  thdt  described 
on  the  N sketch,  and  thereafter  into  that  of  the  SB  sketch. 

Now  certain  inkling  of  stability  becomes  felt.  If  the  spin  is 
sufficiently  great,  the  angular  velocity/!/  of  quasi-nutation  is  also 
great,  and  that  of  quasi-precession  (essentially  the  velocity  M/L  of 
L divided  by  the  magnitude  of  the  qua3i-precession) , small.  The  plane 
of  jm  then  begins  to  catch  up  with  the  quasi-precession.  However, 
since  the  yaw  is  still  growing,  the  velocity  M/L  is  still  increasing. 

It  is  still  possible,  therefore,  that  the  yaw  might  not  catch  up  with 
the  quasi-precession  (via.,  might  not  catch  up  with  L).  It  makes 

only  a theoretical  difference  if  it  does  catch  up  eventually,  but 
only  so  late  that  the  yaw  builds  up  to  a large  value.  If  it  fails 
to  30  catchup  (and  In  practice,  if  it  catches  up  too  late),  the  shell 
is  unstable.  Note  that  throughout  this  NE  regime  the  magnitudes  of 
the  quasi-nutation  and  of  the  quasi-precession  continue  to  grow. 

Generally,  however,  catches  up  with  L,  and  the  situation 

passes  into  that  of  the  E sketch,  i.e.,  an  instantaneous  right-handed 
circular  yawing.  The  quasi-nutation  and  the  quasi-precession  now 
cease  their  growth.  It  is  now  that  the  stability  of  the  shell  is 
finally  and  truly  manifested}  in  fact,  it  now  becomes  of  interest  to 
formulate  the  quantitative  relations  underlying  the  stability  (this 
we  shall  presently  proceed  to  do)* 

/ 

Thereafter,  of  course,  the  situation  passes  into  that  represented 
by  the  SB  sketch}  the  yaw,  the  quasi-nutation  and  the  quasi-precession 
begin  shrinking  in  magnitude . This  continues  through  the  regime 
described  by  the  S sketch  and  into  that  of  the  5W  sketch}  after  which 
L overtakes  again,  and  the  situation  passes  into  the  W sketch  again. 

Throughout  the  SE,  S and  Sd  regime  the  kinetic  energy  due  to  the  cross- 
spin  T|  is  being  converted  into  the  potential  energy  of  yaw,  as  the  shell 
is  brought,  against  the  overturning  torque,  to  face  mare  nearly  into 
the  wind.  Thereafter  the  potential  energy  slips  back  (in  the  1W,  N and 
NB  regimes)  into  the  kinetic  energy  of  cross-spin* 

At  this  stage  we  do  not  care  (as  yet)  to  assert  that  the  shell 
returns  to  its  original  angles  of  yaw  and  that  the  motion  has  a periodic 
character  - even  though  this  is  indeed  the  case  once  our  assumptions 
are  properly  crystallized.  In  reality  there  exist  a number  of  complicating 
factors.  Thus,  the  successive  loops  might  become  smaller,  and  eventually 
become  cusps}  the  epicycle  then  becomes  a ’’rosette*,  then  a wavy  circle, 
and  finally  passes  into  a precession  similar  to  that  described  above* 

Or,  the  loops  might  grew}  the  epicycle  would  then  puss  through  T^} 

thereafter  the  loops  would  encircle  Tp  become  more  nearly  like  coincident 
circles,  and  the  epicycle  passes  into  nutation.  Again,  the  epicycle 
might  shrink  (or  grow)  retaining  its  shape.  What  matters  for  our  present 
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purposes  is  that  this,  roughly  described,  mechanism  of  the  motion  survives 
in  any  case;  and  that  any  changes  in  the  motion  of  A^  can  be  accounted  by 
the  changes  in  the  motion  of  L and  T^. 

STABILITY 


Concept 

A theory  of  yawing  motion  of  a shell  is  essentially  the  theoiy  of 
Stability  of  a shell.  While  much  has  been  written  on  the  concept  of 
stability,  the  situation  still  seems  unsatisfactory  in  spots.  The  word 
conveys  an  idea  that  is  necessary  and  useful,  tut  vague  and  having  many 
ramifications.  The  tradition  of  mathematical  approach  demands  that  the 
word  be  associated  with  a definite  quantitative  criteribn:  for  instance, 
the  shell  is  considered  (in  the  linear  theory)  stable  (in  a certain  sense) 
if  the  so-called" stability  factor"  is  greater  than  1 (and  .then  again,  it 
can  be  less  than  0 - and  also,  slightly  less  than  1 for  spin-stabilised 
rochets).  The  trouble  is,  there  are  so  many  ramifications.  We  have, 
for  instance,  the  static  stability  (fin-stabilisation),  gyroscopic 
stability  (spin-stabilisation),  dynamic  stability  (the  requirement  that 
both  precession  and  nutation  are  damping  out  - distinctly  a concept  of 
linear  theory);  many  ingenious  combinations  of  these  requirements  are  in 
use.  The  subject  of  the  trailing  of  the  shell,  obviously,  is  an  aspect 
of  stability.  Accustomed  to  thinking  in  terms  of  the  linear  theory,  we 
are  wont  to  take  it  for  granted  that  we  can  separate  the  question  of 
stability  from  the  question  of  the  magnitude  of  yaw  (that  is  in  effect 
dismissed  as  a question  of  the  constants  of  integration);  but  it  is 
awkward  to  speak  of  the  stability  of  a shell  which  is  performing  a 
steady  nutation  of  85>  amplitude,  and  it  is  unfair  to  speak  of  the 
instability  of  a shell  which  in  its  whole  trajectory  changes  its  yaw 
from  1 to  2 • In  the  non-linear  theory  the  concept  of  stability  be- 
comes dependent  upon  (±;e.,  inseparable  from)  the  magnitude  of  yaw 
and  other  initial  conditions.  In  the  more  advanced  theories  there 
exists  a need  for  still  other  stability  criteria,  that  have  sot  as 
yet  been  formulated.  Thus,  in  fact,  the  complexity  of  the  problem 
defeats  the  good  intention  of  making  an  intuitive  idea  concrete. 

We  feel,  nevertheless,  that  the  intuitive  idea  of  stability  of 
a shell  is  prior  to,  and  independent  of,  all  quantitative  criteria. 

We  therefore  propose  that  the  essence  of  the  concept  of  stability  la 
a predictions  an  assurance  that  the  yaw  will  not  reach  large  values. 

Naturally,  a prediction  can  be  only  as  good  as  the  information 
on  which  it  is  based.  If  - like  a Maxwell  demon  - we  knew  literally 
all  there  is  to  know  about  the  shell  at  a given  instant  we  would  surely 
be  able  to  predict  exactly  the  behaviour  of  the  shell  for  the  remainder 
of  its  flight.  The  more  advanced  theories  of  exterior  ballistics  can 
do  no  more  than  attempt  to  approach  this  situation.  The  linear  theory 
obviates  the  difficulty  (of  our  ignorance)  by  tha  simple  expedient  of 
resorting  to  clearly  specified  assumptions;  it  is  precisely  this 
apparent  ambition  of  the  linear  theory  which  gives  to  the  student  of 


ballistics  the  false  hope  that  a complete  prediction  is  possible  (and  which 
leads  to  the  natural  - and  at  this  time  apparently  well-entrenched  - demand 
that  any  theory  give  a complete  prediction).  The  fact  is,  the  accuracy, 
and  the  surety,  of  a prediction  depends  upon  the  range  (say,  the  interval 
of  time)  to  which  this  prediction  is  to  apply.  In  most  practical  cases, 
probably,  the  linear  theory  can  give  an  excellent  prediction  for,  say, 
several  periods  of  yaw  (e.g.,  a linear-theory  epicycle  can  be  fitted  to 
well-nigh  every  stretch  of  several  periods  of  yaw).  For  a time  interval 
of,  say,  one  eighth  of  the  period  of  quasi-natation  the  reader  can  make 
a fair  prediction  simply  by  visualizing  one  of  outr  sketches.  But  in  the 
very  general  case  (say,  a highly  complicated  force  system)  a good  prediction 
can  be  made  only  for  a very  short  time  interval. 

What,  then,  is  the  meaning  of  stability?  Obviously,  it  must  be  viewed 
in  relation  to  the  length  of  time  to  which  our  prediction  will  apply.  We 
may  speak  of  “local"  (instantaneous)  stability;  of,  say,  "short-range" 
stability  (the  assurance  that,  if  the  yaw  is  increasing,  the  increase  will 
come  to  a stop  and  will  be  followed  by  a decrease);  and  finally,  of  various 
degrees  of  "long-range"  stability  (an  assurance  that  the  subsequent  maxima 
of  yaw  - second,  third,  perhaps  hundredth, will  not  exceed  a specified  number. 

In  the  regimes  of  SE,  S and  SW  sketches  the  local  stability  is  inherent 
in  the  initial  conditions:  the  yaw  is  decreasing.  In  the  regime  of  HW  sketch 
the  stability,  even  if  it  does  exist  in  some  sense,  is  nc^  manifest.  Even 
if  the  IW  regime  passes  into  the  HE  regime,  the  stability  is  not  certain; 
stability  can  be  implied  only  by  the  assurance  that  the  motion  will  indeed 
reach  the  situation  described  by  the  E sketch,  when  the  point  A^  rolls,  so 

to  say,  over  T^.  Let  us  new  inspect,  quantitatively,  the  situation  at  such 
an  instant. 


Formulation 


Let,  on  the  E sketch,  the  magnitude  of  yaw  (A^T^)  be  5,  the  magnitude 

of  quasi-nutation  (IA^)  be  0,  let  M be  positive  (spin-stabilized  shell)  and 

let  be  stationary  (rectilinear  flight).  The  velocity  r\  of  isflsin^, 

or  (L/B)si n0  (cf.  equation  (2));  the  angular  velocity  of  A^  about  is 

(L/B)sin0/sin6.  The  velocity  of  L is*  M/L;  the  angular  velocity  of  L 
about  is  (M/L)/sin(6  - 0).  The  stability  requirement  is 

( L/B) s in0/s in.5  ^ (M/L)/sin(S  *0), 


and  (since  in  this  example  M is  positive,  and  no  angle  exceeds  n),  can 
be  written  as 


(L/B)sin0/sin5 
(M/L)/sim(6  - 0) 


L sin0  sin(6  - 0)  ^ ' 

B g(5)  staS  » lj 


(3) 


where  M is  written  in  the  form  which  reminds  us  that  it  is  the  function  of 
of  6.  This  is  the  stability  factor  relation  in  a generalized  farm  (although 
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Generally,  it  is  Mcos(L,B)/L  (cf.  Fig.  1*);  but  for  the  E sketch  the  angle 
(L,B)  is  zero. 


the  left-hand side  of  this  inequality  is  not  necessarily  the  stability 
factor,  as  shall  be  explained  presently). 

Linear  Skeleton 


For  the  sake  of  both  simplicity  and  check,  let  us  revert  briefly  to 
the  assumptions  of  the  linear  theory.  ThenL  ^ (because  of  the  smallness 
of  0 - cf.  Fig.  l),  M c p6  (ji  a constant),  and  sines  are  equal  to  their 
angles;  (3)  becomes 


a2  2 
A 

-sir 


(0/6)  (1  - 0/6)  s*  1 


(30 


The  product  (0/6)  (l  - 0/6)  has  a maximum  value  of  l/kj  occurring  at 
0 » 6/2 1 even  with  this  0 , there  remains  to  be  satisfied  the  inequality 


(3») 


which  is  the  definition  of  the  stability  factor  s.  The  physical  significance 
of  the  stability  factor  is  thus  particularly  simples  for  a given  combination 
of  A,  B,  and  n,  the  stability  factor  is  the  maximum  possible  ratio  of  the 

angular  velocity  of  the  plane  of  yaw  to  the  angular  velocity  Of  quasi-pre- 
cession.  This  ratio  occurs  only  at  very  particular  cir cumstances  ; at  the 
maximum  of  yaw,  and  when  the  vector  of  angular  momentum  bisects  the  angle 
of  yaw.  It  is  at  these  circumstances  that  the  stability  is  most  clearly 
manifested.  It  is  interesting  to  note  that  the  epicycle  of  the  linear 
theory  in  these  circumstances  passes  through  the  origin;  i.e.,  this  is 
the  familiar  case  resulting  from  the  initial  conditions  6-0,  d6/dt  / 0. 

Given  s>l,  (3*)  can  be  satisfied  not  only  at  0 • 6/2,  but  also 
in  a certain  neighboring  range  of  0.  At  the  extremes  of  this  range  (3*) 
is  barely  satisfied,  i.e.,  becomes  the  equality 

Its (0/6)  (1  - 0/6)  « 1, 
which  can  be  rewritten  as 

(0/6)2  - (0/6)  + lAs  - 0,  (3A) 


and  has  the  two  solutions, 

(0/6)l,2  * \ (1  i\/l  “ Vs)  - (1  1 *)/2,  say,  (3B) 


These  two  extremes  are,  of  course,  nutation  and  jrecession.  In 
nutation  L is  closer  to  T^,  or  0/6  is  more  than  1/2;  A^  moves  fast 

because  of  the  large  0,  and  L just  barely  keeps  up  with  A^.  In 

precession  L is  closer  to  A^j  and,  be  cause  of  the  small  0j  A^  moves 
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only  slowly,  and  just  barely  keeps  up  with  L.  The  instantaneous  motion 
of  the  body  in  both  cases  is  rotation  about  co  at  the  same  angular  velocity 
a>;  or,  the  motion  of  in  both  cases  is  rotation  about  L at  the  same 

angular  velocity  Pi*  But  the  reader  should  remember  that  these  angular 
velocities  describe  the  rotation  (of  the  body,  or  of  the  8 vector”  of  quasi- 
nutation)  about  the  instantaneous  (i*e«,  stationary,  as  it  were)  position 
of  oo  and  Lr  and  that  the  points  co  and  L are  moving.  The  situation  i3 
simple  enough  with  quasi-process  ion  (L  moves  with  velocity  M/L),  or  with 
the  rotation  of  the  body  (it  is  rolling),  but  with  regard  to  the  quasi- 
nutation  the  situation  is  not  quite  the  same  sort  of  thing  for  which  one 
is  accustomed*  The  concept  of  the  angular  velocity  of  the  quasi-nut  at ion 
could  have  been  Introduced  as  the  angle  between  two  subsequent  positions 
of  quasi-nutation  divided  by  the  time  increment ;*  this  was  not  done 
because  such  an  angular  velocity  would  not  be  the  same  as  the  constant  fl , 
and  would  not  be  as  convenient;  different  such  velocities  can  be  produced 
by  the  saae*A. 

Tie  variation  of  the  shape  of  epicycle  in  the  range  (l  - 6)/Z 
< 0/6<(l  + 6)/Z  is  very  interesting.  If  we  decrease  0/6  below  l/2t  the 
loops  begin  to  miss  the  center,  become  smaller,  and  epicycle  passes  through 
a rosette  to  wavy  circle  and  into  precession;  if  we  increase  0/6  ever  1/2, 
the  loops  start  going  around  the  origin,  increase,  and  the  epicycle  passes 
into  nutation.  In  either  case,  the  variation  of  the  magnitude  of  yaw 
diminishes,  until  it  finally  disappears  in  the  corresponding  circular  motion. 
Thus  the  rapid  variation  of  the  magnitude  of  yaw  (as  in  the  epicycle 
passing  through  the  origin,  with  0/6  =*  1/2)  is  a token,  or  a manifestation, 
of  stability. 

Outside  of  this  range  of  0/6  the  inequality  (3*)  is  not  satisfied, 
and  the  stability,  even  if  it  does  exist  in  some  sense,  Is  not  manifest 
in  the  regime  of  the  E sketch*  This  sketch  then  represents  a minimum, 
rather  than  a maximum,  of  yaw.  With  the  smaller  0/6,  the  epicycle  will 
be  a wavy  circle;  with  the  larger  0/6,  the  resultant  epicycle  will  loop 
over  the  origin* 

A very  important  borderline  case  is  that  of  s * 1*  Precession 
and  nutation  then  merge,  and  only  a bare  manifestation,  or  a minimum,  of 
stability  (viz*,  circular  motion)  becomes  possible;  and  that,  only  at 


In  fact,  In  the  linear  theory,  when  nutation  is  considered  a vectorial 
component  of  yaw,  this  is  exactly  how  the  angular  velocity  of  nutation 
is  defined;  it  is  useful  there,  for  such  angular  velocity  then  happens 
to  be  constant*  Presently,  when  we  shall,  relax  the  assumption  of  the 
stationary  T^,  we  shall  meet  the  same  difficulty  with  the  quasi- 
precession* 
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the  most  advantageous  condition,  0/6  c l/2  . Note  that  such  motion  is 
still  possible  at  all  (small)  angles  of  yaw. 

^ore  General  Gases 


With  this  understanding  of  the  skeleton  of  the  "short-range®  (in 
this  case,  gyroscopic)  stability  we  may  relax  our  assumptions  of  linearity, 
and  return  briefly  to  the  more  general  problem,  i.e.,  to  (3)*  New  we 
have  more  exactly  (subject,  of  course,  to  other  assumptions  here  implied) 

I = Acol./cos0  (cf . Fig.  l),  and  M **  |isin6  (where  p.  so  defined  is  not 
necessarily  a constant,  but  is  certainly  more  nearly  so  than  M).  Then 
(3)  becomes 


2 2 

A “1  sin/  sin (6-0)  ^ , 

5IC5T  — L 21  ^ 

^ ' cos  0 son  6 


( k ) 


which  we  shall  write  as  sf(0,6)^l,  with  a x A^aa^A.B^i.(6)  and  f e 

iisim0sin(6-^O/cos^0sin^6.  While  (1*)  is  still  similar  to  its  linear 
analog  (3')>  the  important  difference  is  that  with  the  larger  angles  the 
trigonometric  part  of  (k),  cr  f/k,  can  have  values  greater  than  l/k>  and 
the  stability  requirement  is  thereby  made  less  strigBnfc. 


The  very  interesting  function  f(0,6)  is  given,  as  a contour  plot  in 
the  (0,6)  plane  for  O<0<6<n/2,  on  Figure  6.  The  most  important  region  of 
this  surface  is  the  region  of  small  0 and  6,  where  the  two  contour  lines  of 
f ■ 1 merge;  this  is  the  domain  of  the  linear  theory.  In  this  region  the 
function  has  a practically  horizontal  crest,  or  watershed;  this  crest 
extends,  as  a razor  wedge,  right  to  the  point  0 6 ■ 0,  which  represents 

a steady  nose-on  flight.  However,  on  the  boundary  lines  radiating  from 
the  crigin  the  function  is  zero,  and  the  manifestation  of  stability  is 
impossible.  In  particular,  the  abscissa  (0/6  « 0)  represents  the 
situation  when  the  axis  of  a spinning  shell  is  momentarily  stationary 
(as  inQthe  .cusp  of  a rosette),  and  1 is  sweeping  across  the  plane  of  yaw. 

The  h$  line,  0/6  *=  1,  represents  the  situation  when  I sweeps  across 
the  trajectory,  and  therefore  rotates  at  an  infinite  angular  velocity, 
as  it  were,  about  the  trajectory.  In  the  linear  region  the  crest  lies 
along  the  line  0/6  = 1/2;  on  the  extension  of  this  line  (which,  for 
moderate  yaws,  is  not  far  from  the  true  crest)  the  function  has  the 


The  mathematical  linear  theory  shows  that  when  two  independent  solutions 
merge  into  a single  one,  as  here,  there  arises,  as  the  second  possible 
solution,  the  "secular"  solution;  in  this  case,  it  is  an  expanding 
Archimedean  spiral.  When  s becomes  less  than  1,  both  the  circular  and 
the  Archimedean  solutions  pass  into  logarithmic  spirals  (respectively, 
shrinking  and  expanding).  The  blend  between  various  types  of  solutions 
can  be  inspected  as  a matter  of  initial  conditions  (constants  of 
integration).  For  instance,  with  arbitrary  initial  conditions  and  with  s 
barely  over  1.,  the  NE  regime  will  be  veiy  long,  and  the  maximum  angle 
of  yaw  very  iGfrge;  with  s of  1 or  less  the  NE  regime  will  continue  until 
the  yaw  reaches  values  so  great  that  the  linear  theory  no  longer  applies. 
We  have  avoided  such  difficulties  by  the  simple  device  of  considering 
only  the  initial  conditions  of  the  E sketch. 
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value  l/cos^(6/2),  which  (for  these  yaws)  is  only  slightly  smaller  than 
l/cos6. 

The  point  0 = 6*=  n/2  represents  a pure  cartwheeling  of  the  shell, 
i.e. , cartwheeling  without  any  axial  spin.  At  this  point  the  function 
f may  have  any  value  from  zero  (since  it  is  on  the  line  0/6  = l)to 
infinity  (since  it  is  also  on  the  line  6 = n/2,  which  represents  the 
cartwheeling  of  the  axis  of  the  shell,  and  on  this  line  f becomes  simply 
iitan0).  Indeed,  at  such  conditions  manifestation  of  stability  rather 
loses  its  meaning. 

It  will  be  seen  that  the  true  crest  of  our  surface  approaches  this 
"infinite  peak"  (the  singularity)  in  a rather  devious  way,  first  going 
outside  of  the  limits  of  the  triangle  0 < 0 < 6 < n/2 , It  is  important  to 
note,  however,  that  this  "peak18  can  be  reached  without  going  outside  of 
this  triangle  ; i.e.,  all  values  of  f are  represented  in  this  triangle. 

This  means  that  for  any  shell  which  is  unstable  in  the  linear-theory 
sense  there  may  exist  combinations  of  0 and  6 such  that  this  shell  will 
still  manifest  stability. 

The  stability  condition  can  then  be  visualized  as  fellows.  For  a 
given  s (which  we  may  new  term  "linear  stability  factor" ) we  may  Imagine 
on  Figure  6 a point  with  map  coordinates  0 and  6 and  at  the  height  l/s  • 

If  this  point  is  above  the  f-surface,  the  stability  condition  (it)  can  not 
be  satisfied  (the  regime  of  the  £ sketch,  even  if  does  exist  instantaneously, 
can  not  continue,  and  will  pass  into  the  regime  of  the  ME  sketch).  If  this 
point  is  in  the  f-surface  (on  the  contour  line  f ■ 1/s),  the  circular  yawing 
will  continue  (the  stability  is  barely  manifested)  • If  this  point  is  under 
the  f-surface  (inside  the  loop  of  the  contour  line  f ■ l/s),  the  stability 
is  clearly  manifested  (the  regime  of  the  E sketch  will  pass  into  that  of 
the  SE  sketch,  and  the  yaw  will  presently  start  decreasing). 

Naturally,  the  possible  circular  yawing  regimes  might  be  divided 
(say,  by  the  crest  of  the  function  f)  into  those  of  nutation  and  precession; 
the  rates  of  nutation  and  precession  can  thus  be  computed,  either  as  the 
rate  of  the  plane  of  yaw  (L/B)sin0/sin6,  or  as  the  rate  of  the  quasi- 
precession  (M/L)/sin(6  - 0).  We  may  thereby  relax  new  the  requirement 
(mentioned  above  as  a matter  of  convenience  of  the  presentation)  that  M 
be  small.  These  rates,  however,  are  no  longer  independent  of  the  initial 
conditions  (the  magnitude  and  the  angular  rate  of  yaw).  Moreover,  these 
two  regimes  are  no  longer  separated^ on  Figure  6,  by  the  regime  of  the 
circular  yawing  motion  of  a barely^projectile  (as  they  are  in  the  linear 
theory);  rather,  nutation  and  precession  can  now  be  blended  across  the 
crest  of  the  f-surface.  The  multiplicity  of  possible  solutions  in  the 
region  of  this  blend  corresponds  to  the  fact  that  in  the  linear  theory 
a barely  stable  shell  may  yaw  circularly  at  all  (small)  angles  of  yaw. 

It  may  be  noted  on  Figure  6 that  the  true  crest  of  the  f-surface 
does  not  represent  exactly  the  minimum  values  of  6 at  which  a given 
linearity-unstable  shell  may  manifest  its  stability.  The  interesting 
question  of  what  these  values  are  has  a particularly  simple  answer: 
the  locus  of  the  minima  of  6 can  be  found  by  finding  the  maxima  of  f 
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ia  the  planes  of  6 = constant.  The  0-dependent  part  of  (U)  can  he  written 
as 

2 

sin0(sin5cos0  - cos6sin0)/cos' 0 =»  cos6  tan/  (tanb  - tan/), 

which  is  of  the  form  ax(b  - x),  and  has  a maximum  at  x = b/2,  i.e.,  at 
tan/  « (tanS)/2.  Thus,  the  vector  of  angolas  momentum  then  bisects  the 
tangent  of  the  angle  of  yaw  erected  on  the  axis  of  the  3heH  (Figure  7a), 
rather  than  the  angle,  as  in  linear  theory*  This  locus  merges,  of  course, 
with  the  line  0/5  *»  1/2*  At  these  conditions  the  value  of  f is  3imply 
l/co36.  Therefore,  the  generalization  of  the  stability  condition  (3“) 
of  the  linear  theory  is 


s/cos6>l,  (U1) 

with  the  proviso  that  at  the  critical  6 stability  nay  be  manifested 
only  if  the  yawing  is  instantaneously  circular  and  ft  has  a particular 
value*  This  proviso  corresponds  to  the  fact  that  in  the  linear  theory 
a barely  stable  shell  may  manifest  stability  only  at  particular  initial 
conditions,  viz*,  at  0/5  ■ 1/2. 

In  fact,  in  advanced  texts  on  ballistics  the  stability  factor  is 
defined  with  cos6  in  the  denominator,  i.e*,  as  the  left-hand  side  of 
(U*)  rather  than  by  (3H ) j we  might  then  speak  of  the  H non-linear" 
stability  factor.  Ills  is  quite  a different  thing  from  our  linear  stability 
factor,  and  no  longer  can  be  viewed  as  a definite  property  of  the  shell, 
independent  of  the  initial  conditions  (as  the  linear  stability  factor 
is  usually  considered)  3 rather,  it  should  he  viewed  as  the  maximum  possible 
ratio  of  the  angular  velocities  of  yaw  and  quasi-precession. 

The  requirement  of  a sufficiently  high  yaw  (6)  may  be  interpreted  as 
a requirement* of  a sufficiently  large  quasi-nutation  (0)  and  a sufficiently 
large  quasi-precession  (6  - /).  It  might  complete  a certain  symmetry  of 
our  arguments  if  we  inquire,  briefly,  if  loci  of  such  minima  exist,  la 
fact,  the  locus  of  minima  of  0 is  quitg  pronounced  on  Figure  6j  it  exists 
from  the  linear  domain  to  6 of  some  56°.  Its  equation  is  tan(6  - 0)  - (tanS)/2 
i.e.,  for  stability  0 must  be  sufficiently  large  at  least  to  bisect  the 
tangent  of  the  angle  of  yaw  erected  on  the  trajectory  (Figure  7b)*  For 
larger  5 this  curve  represents  the  maxima,  instead  of  the  minima,  of  0 on 
the  contour  lines  of  f . The  locus  of  the  minima  of  6 - 0 is  bagely 
noticeable  in  the  region  from  the  linear  domain  to  6 of  some  25°,  and  seems 
to  be  of  no  particular  interest. 

The  usefulness  of  the  stability  factor  lies  in  its  relation  to  the 
nutational  and  precessional  rates.  In  the  linear  theory  these  rates  are  - 
from  (2),  Figure  1 and  (3B)  - simply 

^ /5  <*fl0/5  m ( Ja^/B)  (0/6)  m (Attj/; 2B)(l  W 1 - l/a ) (5) 


* For  any  s < 1. 
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For  any  circular  motion  the  nan-linear  analog  of  (3A)  can  be  shown  to 
be 

o 

tan  0 cctS  - tan0  + sinS/lta  ■=  0,  (ItA) 

so  that  the  analog  of  (3B)  is 

tan0/taa6  » (l  + Vl  - coa6/s)/2,  (liB) 

and  the  circular  rates  are 

"!}  /ainfi  =?fl  ain0/sin6  = (A<o^/B)(tan0/sin6)  = (iW^/2B)  (l  + ^1-c036/s)/cos6*  (6) 
Example 

The  foregoing  discussion  was  mainly  the  recital  of  the  more  interesting 
conclusions  from  the  theory  of  spinning  top  - for  which  case  the  assumption 
of  a constant  (t  is  particularly  justified*  While  the  analogy  between  the 
shell  and  the  top  is  quite  traditional,  a comparison  between  the  linear 
mathematical  analog  and  the  every-day  experience  with  a top  is  met  rather 
less  frequently;  such  a comparison  might  be  rather  illuminative* 

In  the  actual  top  the  nutation  damps  out  extremely  rapidly  (by 
processes  of  which  more  anon,  in  connection  with  Magnus  torque;,  so  that 
in  practice  it  is  observed  seldom  if  ever*  The  precession  also  damps 
oat  quite  rapidly;  both  the  damping  and  the  angular  rates  of  precession 
are  readily  sensed  when  we  see  a hard -spinning  top  right  itself  and  "go 
to  sleep" . let  the  reader  will  recollect  that  the  top  often  fails  to 
regch  the  sleeping  position,  and  continues  a rather  steady  "precession"; 
also  as  the  spin  of  a sleeping  top  is  dying  out,  there  invariably  develops 
the  same,  relatively  slow,  precession  - which  speeds  up  as  the  yaw 
increases,  just  before  the  top  falls*  These  familiar  "steady"  motions, 
apparently,  are  not  the  precession  in  the  customary  ballistic-  sense 
(which  is  a result  of  disturbance*  expressed  by  a solution  of  a 
homogeneous  differential  equation};  rather,  they  are  the  motions  of  a 
line arly-un* table  top,  which  is  barely  stable  (has  the  non-linear 
stability  factor  a/cos6  of  1)  at  the  given  angle  of  yaw*  Any  oscillations 


The  apparently  peculiar  behaviour  of  the  last  expression  in  (6)  at 
6 » n/2  can  be  readily  accounted  for  by  inspecting  the  legitimacy  of 
our  mathematical  manipulations  at  this  6.  The  processional  rate 
(which  appears  indeterminate)  then  is  simply  jWtan0/B,  i*e*,0Q . That 
the  mtational  rate  is  not  infinite  can  be  seen  as  follows^  From 
(UB)  or  Figure  6 it  is  seen  that  natation  at  6 ■ n/2  is  possible  only 
if  0 ■>  n/2,  i*e*.  if  the  axial  spin  to*  is  zero,  and  the  apparently 
Indeterminate  (6)  then  amounts ^to  7]  again,  - with  the  additional 
provisb  that  if  the  motion  is  to  be  circular,  as  assumed,  M must  be 
zero*  The  trick  difficulties  of  this  kind  are  inspected  best  by 
examining  the  physical  significance  of  the  problem  (e.g.,  Figures  2 
and  3)j  rather  than  the  legitimacy  of  algebra. 
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of  this  angle  rapidly  damp  out,  apparently  by  processes  similar  to 
the  damping  of  the  ballistic  nutation  and  precession;  in  case  the 
non-linear  stability  factor  is  less  than  1,  the  top  merely  builds  up 
its  yaw  till  it  is  non-linear ly  barely  stable  at  same  larger  yaw;  the 
speedup  of  the  precession  (in  spite  of  the  decrease  of  co^)  Just  before 

the  top  falls  is  the  consequence  of  the  presence  of  cos6  in  the  denominator 
of  (6).  Thus  it  is  the  usually-neglected  border-line  case  of  the  linear 
theory,  rather  than  the  theoretical  epicycle,  which  is  the  ancestor  of 
the  commonly-observed  motions  in  the  non-linear  regime.  It  is  extremely 
easy  to  confuse  this  motion  of  non-linear  regime  with  the  ballistic 
(linear)  precession;  i.e.,  to  consider  it,  erroneously,  a matter  of  initial 
conditions.  A shell  with  an  insufficient  linear  stability  factor  may 
not  exactly  "tumble”;  it  is  still  "stable”  in  'h  sense.  It  is  Just  that 
its  yaw  is  undesirably  large. 


PREDICTION  OF  THE  YAWING  MOTION 

In  the  tradition  of  mathematical  approach  the  stability  condition 
stands  somehow  apart  from  the  solution;  although  it  is  understood  that 
it  is  also  a kind  of  a short-cut  summary  of  the  character  of  the  solution. 
We  interpreted  the  stability  condition  as  an  attempt  at  the  prediction 
(i.e.,  the  solution)  itself;  and  we  found  that  a consideration  of  an 
instantaneous  situation  may,  in  a number  of  important  special  cases, 
yield  a modicum  of  a long-range  prediction.  We  avoided  a quantitative 
discussion  of  the  local  stability  (viz.,  of  the  question  whether  a pre- 
diction for  a short  length  of  time  augurs  well  or  ill  for  the  decrease 
of  yaw);  and  discussed  quantitatively  only  the  short-range  stability, 
which  is,  as  it  were,  the  local  stability  at  the  instant  of  the  first 
maximum  of  yaw.  £ut  to  get  to  the  first  maximum  of  yaw  from  arbitrary 
initial  conditions,  as  well  as  to  get  from  the  first  maximum  to  the 
second,  third,  etc.,  is  a different,  and  a more  difficult,  problem; 
in  fact,  this  is  the  problem  which  has  been  solved  analytically  only 
in  the  most  simple  cases  (such  as  the  linear  theory,  or  the  case  of 
the  M-torque  acting  alone),  and  which  requires,  generally,  modern 
high-speed  computing  machinery. 

In  the  linear  theory  the  general  solution  is  obtained  by  the  ' 
superposition  of  the  simple  nutation  and  precession;  and  the  super- 
position arises  naturally,  as  a property  of  solutions  of  a linear 
differential  equation.  As  a matter  of  a mere  exercise,  or  of  an 
inquiry  hew  far  we  can  go  without  calculus,  the  following  interpretation, 
or  a derivation  of  the  superposition  might  be  of  insterest. 

For  any  yaw  5 (which  in  the  linear  theory  is  simply  a complex 
number)  there  exist  two  poesSiLenafces,  a^d  at  which  the 

motion  would  be  the  simple  (circular)  motion;  for  each  one  of  these 
possible  rates  the  given  £ must  be  decomposed  into  quasi-precession 
and  quasi-nutation  in  a definite  way  - which,  as  we  have  seen  by  (3A), 
is  independent  of  £ • The  actual  instantaneous  rate  f of  \ , generally. 
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is  neither  one  of  these.  Now,  two  complex  numbers,  and  Cp,  can  easily 
be  found  such  that  C + C - 1,  and  t » G C + C t Obviously,  t may 
be  decomposed  into  Cq  £ and  £ j a brief  geometrical  or  algebraic  con- 
sideration, detailed  on  Figure  8,  will  show  that  the  quasi-nutation  and 
quasi-precession  may  be  decomposed  similarly  (e.g.,  on  Figure  8, 

QP  = + CpQPp)  j the  vectors  QPn  and  QPp  are  designated  on  Figure  8 

by  L and  L • Also,  the  projected  vector  of  the  torque  M may  be 

r 

decomposed  in  the  same  way. 


Now,  it  does  not  matter  which  component  of  torque  acts  on  which 
component  of  angular  momentum,  as  long  as  the  total  M acts  upon  the 
total  L.  Furthermore  (but  this  holds  only  in  the  linear  theory),  it 
does  not  matter  how  the  quasi-nutation  is  decomposed.  All  of  its 
components  are  rotating  (about  the  corresponding  stationary  positions 
of  L)  at  the  same  angular  velocity  -ft,  and  the  velocity  f (or 7) ) of  A^ 

is,  in  any  case,  the  vectorial  sum  of  the  velocities  produced  by  the 
components  of  quasi-nutation.  We  therefore  can  assign,  arbitrarily, 
the  following  relations.  Let  the  torque  produced  by  the  component 

of  yaw  affect  the  CnQPn  component  of  quasi-precession  QPj  and  let  the 

angular  momentum  component,  represented  by  that  component  of  quasi- 
pbecession,  produce  - by  means  of  the  corresponding  component  CnQNn 

of  quasi-nutation  QN  - the  component  Cn  £ n of  j with  similar  relations 

holding,  of  course,  among  the  components  with  the  p-subscript.  We  can 
then  note  that  among  the  components  with  each  subscript  the  motion  is 
known,  i.e.,  predictable,  circular  motion.  Therefore  the  G's  will 
remain  constant,  and  the  general  solution  is  a superposition  of  precession 
and  nutation. 


While  this  intrepretation  is  particularly  obvious  in  a highly 
simplified  case  (rectilinear  flight,  constant  velocity,  only  M present), 
it  is  not  vitiated  if  we  pass  to  the  general  (linear)  case,  fle  need 
only  to  re-label  the  components  Cn£  and  Cp  £as  Cn£n(t)  ^ C f (t), 

introduce  the  yaw  of  repose  and  the  more  complicated  system  of  torques 
and  forces,  (all  of  which,  however,  must  be  linear  homogeneous  functions 
of  (*  and  £ ),  introduce  the  motion  of  T^,  etc.,  and  re-trace  the 

argument  more  slowly  and  more  laboriously.  In  fact,  the  component 
motions  £ Q(t)  and  £^(t)  may  be  quite  complicated:  e.g.,  such  parameters 

of  motion  as  axial  spin,  velocity  and  stability  factor  may  all  change  - 
provided  they  change  in  the  same  way  for  both  £ and  £ 


The  latter  requirement  rules  out,  in  the  linear  theory,  such  phenomena  as 
instability  of  spin  (BEL  668),  cases  when  the  yaw-drag  affects  the  trajectory, 
etc.;  those  are  domain  of  non-linear  theories.  It  is  not  quite  idle  to  remark 
that  the  component  motions  need  not  even  be  precession  and  nutation  l They  can, 
equally  well,  be  any  two  linearly  independent  epicycles.  Such  procedures  are 
indeed  often  resorted  to,  when  various  linear  perturbations  of  flight  data 
must  be  carried  out. 
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Once  the  superposition  is  accepted,  such  facts  as  the  absence  of 
damping  with  M alone  acting  on  the  shell  (and  with  M being  a function 
of  6 alone),  as  well  as  the  complete  prediction  of  motion,  follow 
readily.  In  this  connection  it  is  well  to  note  the  distinction  between 
the  dynamic  and  the  long-range  stability.  Consider  a shell  whose 
nutation  damps  rapidly,  but  whose  precession  grows,  albeit  very  slowly. 
In  the  clean-cut  tradition  of  the  mathematical  approach,  such  shell  is 
automatically  classified  as  dynamically  unstable.  Our  approach  seems 
(in  that  case)  more  tolerant,  or  more  nearly  attuned  to  the  reality. 

If  the  initial  yawing  motion  is  mostly  nutation,  we  would  say  this 
shell  has  a long-range  stability  for  so-mamy  periods  of  yaw,  or  that 
the  instability  would  not  be  manifested  until  such-and-such  range,  or, 
again,  that  the  stability  depends  upon  the  manner  of  launching.  This 
is  certainly  no  different  from  what  would  usually  be  done:  it  is  just 
that  our  "long-range"  stability  corresponds  to  the  final  calculations 
of  the  linear  theory  - or,  to  proof -firing  results  - rather  than  to 
the  abstraction  of  dynamic  stability. 

An  attempt  to  systematize  the  resultant  possible  multiplicity  of 
epicycles  is  given  on  Figure  17.  It  is  amusing  (and  from  our  "non- 
mathematical"  point  of  view,  chastising)  to  observe  that  in  practice 
the  simplest  way  of  drawing  these  curves  has  been  to  construct  the 
appropriate  differential  equations,  to  compute  the  appropriate  initial 
conditions , and  to  solve  these  equations  numerically  on  an  electronic 
differential  analyzer  (GEDA)  , 

Even  allowing  that  the  general  non-linear  case  is  not  solvable 
analytically,  it  is  of  interest  to  inquire  whether  there  might  not  be 
some  simple  non-linear  cases,  e.g.,  in  the  case  of  H acting  alone,  and 
being  a function  of  5 alone.  Unfortunately,  the  theory  furnishes  no 
simple  and  useful  answers  there.  Even  if  we  are  to  prove  the  very 
relevant  fact,  that  the  yawing  motion  in  such  a case  is  not  damped,  it 
is  difficult  to  keep  up  our  pretence  of  avoiding  the  calculus.  * 

Let  us  now  inspect  the  consequence  of  the  various  components  of  the 
force  system  acting  upon  the  shell  - qualitatively  in  the  general  case, 
and  quantitatively  in  the  linear  case. 


This  work  was  done  by  Mr.  E.  E.  Bomber ger . 

A a was  mentioned  before,  the  general  solution  has  been  worked  out, 
but  is  cumbersome.  In  the  case  of  spinning  top  (constant  n.)  the 
solution  is  in  terms  of  elliptic  integrals. 

Yet  we  might  attempt  to  outline  the  proof,  very  roughly,  as  follows. 
Since  M is  perpendicular  to  the  shell  axis,  the  axial  spin  cannot  be 
affected;  and  since  M is  perpendicular  to  trajectory,  tile  component  of 
L along  trajectory  cannot  be  affected.  It  can  belfiea  shown  that  since 
M depends  upon  6 alone,  there  exists  a potential  energy  of  yaw,  and 
therefore  the  system  is  what  is  known  in  mechanics  as  conservative; 
and  such  systems,  generally,  possess  constant  total  energy  and  represent 
oscillations  "Tn  a potential  well",  with  constant  amplitudes.  Such 
vague  "reasoning"  is  not  without  its  legitimate  .place  in  ballistics; 
in  fact,  it  is  more  logical  than, for  instance,  an  insistence  on  a 
large  stability  margin  of  a mortar  shell!  This  intuitive  groping  for 
"constants  of  motion"  can  be  recognized  as  a primitive  form  of  the 
Hamiltonian  mechanics. 


FORCE  SYSTEMS  CF  BALLISTICS 


la  recent  years  the  fact  began  to  be  recognized  that,,  the  force  systems 
used  in  exterior  ballistics  were  rather  oversimplifications  of  reality.* 
Nevertheless,  for  practical  purposes  an  idealization  is,  of  course,  in- 
dispensable; it  is  just  that  a realization  of  just  what  has  been  assumed, 
or  taken  for  granted,  must  new  be  kept  in  mind  more  prominently. 

The  first  restrictive  assumption  is  that  of  Synge’s  aerodynamic 
hypothesis:  that  the  force  system  (viz.,  total  force  and  total  torque) 
can  be  specified  by  specifying  the  vectorial  velocity  and  angular  velocity 
of  the  shell?*  The  next  systematization  is  the  Maple -Synge  analysis  of  the 
consequence  of  the  symmetry  of  the  shell.’  This  systematization  is 
particularly  useful  in  investigating  the  border-line  between  the  linear 
and  non-linear  theories.  It  restricts  somewhat  the  number  of  possible 
non-linearities  which  ought  to  be  considered.  The  usefulness  of  this 
systematisation  in  the  general  nan-linear  regime  is,  unfortunately, 
limited,  for  two  reasons.  It  is  difficult  to  design  experiments  so 
that  all  probably  relevant  Maple-Synge  coefficients  can  be  measured; 
and  even  if  they  were,  the  mathematical  utilization  of  them  would  be 
extremely  laborious. 

There  are  two  known  ways  in  which  a further  necessary  simplification 
of  Maple-Synge  theory  may  proceed.  The  most  important  one,  of  course, 
is  the  linearization  ***af  the  Maple-Synge  theory.  This  excludes  a number 
of  coefficients  which  are  known  to  be  relevant,  and  includes  same  whose 
relevance  seems  rather  marginal.  The  most  elegant  feature  of  the  linear- 
ization scheme  (in  addition  to  the  fact  that  it  dovetails  with  the  super- 
position of  the  linear  theory)  is  that  it  furnishes  a convenient  coordinate 
system  for  decomposition  of  the  force  system;  forces  and  torqes  dependent 
upon  6 are  decomposed  along  and  across  the  plane  of  yaw,  and  those  dependent 
upon  T)  are  decomposed  along  and  across  the  arrows  of  7)  in  our  sketches.  This 
gives  rise  to  the  well-known  Kelley-McShane  matrices  of  forces  and  torques. 


So  that  a number  of  extant  ballistic  texts  now  appear,  in 'retrospect', 
rather  didactic. 

This  assumption  seems  to  be  rather  peculiar  to  ballistics.  In  the 
advance^  aeronautics  :the  attempts  to  rescind  this  assumption  dates 
to  the  1920's.  Cf , ERL  $82,  "Dynamic  Measurement  on  8lmm,  M 56  Shell 
in  NBS  Wind  Tunnel",  for  reference  and  far  brief  discussion. 

Chronologically,  of  course,  the  linearized  ( Niels en-Synge)  version 
of  Maple-Synge  theory  proceeded  the  more  complete  theory. 
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which  will  be  discussed  in  detail  presently.  Its  -usefulness  as  a 
classification  scheme  survives,  to  some  extent,  even  in  the  non-linear 
theory. 

The  second  known  way  is  the  adoption  of  circular  yawing  motion 
(rather  than  small  yaw)  as  the  basic  form  of  yawing  motion.  The 
most  troublesome  Maple -Synge  coefficients  are  then  excluded,  and 
experiments  can  be  readily  designed  by  which  all  terms  of  the  Kelley- 
McShane  matrix  can  be  measured  - whereby  this  matrix  can  be  considered 
non-linear,  and  merely  an  aid  in  classification.  The  mathematical 
utilisation,  unfortunately,  remains  difficult. 

What  is  probably  wanted  first  is  an  application  of  the  Maple-Synga 
(or  same  other)  theory  of  aerodynamic  forces  acting  on  a given  surface 
of  the  shell  (and  independent  of  the  mass  distribution  within  the  shell) 
to  a class  of  "possible",  or  at  least  "probable®,  yawing  motions  - if 
such  classes  can  be  defined.  Circular  yawing  motion,  for  instance,  is 
useful  because  it  is  a sort  of  approach  to  such  a class.  In  fact,  this 
application  is  achieved,  perforce,  in  free-flight  experiments.  It  should 
be  noted,  however,  that  the  determination  of  the  aerodynamic  forces  from 
the  free-flight  data  is  not  a simple  process;  generally,  it  depend* 
heavily  upon  an  existence  of  a sufficiently  idealised  mathematical  theory. 
Such  a theory  must  start  from  a certain  assumed  force  system,  or,  as  it 
were,  from  a wind-tunnel  point  of  view.  It  is  hoped  that  this  gap,  as  it 
were,  between  the  possible  free-flight  technique  and  the  complete  non- 
linear Maple-Synge  theory  will  be  eventually  bridged  by  extensive  use 
of  high-speed  computing  machinery,  and  by  elaborate  wind-tunnel 
instrumentation  which  would  reproduce  circular, elliptical,  spiral, 
epicyclic,  etc.,  yawing  motion. 

These  vague  generalities  have  cited  here  with  the  object  of  showing 
to  the  student  of  exterior  ballistics  the  present  groping  state  of  the 
art:  it  is  difficult  not  only  to  build  a general  non-linear  theory, 
but  even  to  formulate  a general  non-linear  force  system;  and  this  even 
without  considering  asymmetries  of  the  shell,  liquidity  of  the  filler, 
transitions  from  laminar  flow,  separation  of  the  flow,  intermitt ency 
of  the  turbulent  wake,  etc.  Small  wonder,  then, that  the  linear  theory 
remains  a point  of  diminishing  returns! 


DAMPING  TQBQTE 

Next  to  M in  importance  is  the  component  of  torque  known  as 
"damping"  torque  and  designated  by  H,  Specifically,  this  torque  is 
defined  by  an  understanding  that  indirection  it  opposes  the  angular 
velocity  r]  of  the  shell's  axis,  and  in  magnitude  it  depends  essentially 
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cm  7)  • Ip  the  linear  theory  it  is  assumed  to  be  proportional  to  7\ 

(i.e. , to  £ ),  and  is  therefore  mathematically  analogous  to  viscous 
damping  of  a harmonic  motion,*  The  choice  of  the  term  suggests  that  a 
hope  must  have  been  entertained  that  a torque  so  defined  will  account, 
in  the  main,  for  the  damping  of  yawing  motion.  In  fact,  this  is  very 
nearly  the  case  - although  the  damping  of  the  yawing  motion  involves 
a number  of  phenomena  other  than  H« 

On  Figure  9 H appears  as  an  arrow  directed  from  L to  A- ; i,e,,  it 
tends  to  shrink  the  quasi-nutation  LA^  by  pulling  1 toward  Its 

effect  upon  the  quasi-precession  depends  upon  where  T^  is  with  respect 

to  L and  A^;  and  in  particular,  in  any  circular  yawing  motion  in  which 

X is  between  and  A^,  E tends  to  stretch  the  quasi-precession. 

This  simple  formulation,  however,  has  led  to  one  of  the  pitfalls 
(mentioned  in  the  introduction),  which  stunted  for  a while  the  development 
of  our  "non -mathematical®  approach.  For  instance,  in  the  linear  theory- 
precession  and  nutation  are  substantially  circular  yawing  motions,  at 
a constant  0/5,  i.e.,  with  a constant  proportion  of  quasi-nutation  and 
quasi-precession.  The  theory  shows  that  E damps  the  nutation,  and  causes 
a growth  of  precession.  Thus  in  nutation  H seems  to  cause  a shrinkage 
of  quasi-precession,  and  in  precession  H seems  to  cause  a stretching 
of  -quasi-nutation  - in  a seeming  contradiction  to  our  simple  formulation. 

In  fact,  a stimulating  question  was  once  asked  - In  a circular  yawing 
motion,  how  does  the  torque  E distinguish  between  nutation  and  preoession? 

These  pitfalls  are  sinply  mistakes  that  are  easy  to  fall  into,  but 
also  easy  to  guard  against.  It  is  easy  to  forget  that  a shrinking  or 
stretching  of  quasi-nutation  and  quasi-precession  may  be  accomplished 
by  the  torquQ.H,  as  well  as  by  H.  It  is  also  easy  to  assume  that  a yawing 
motion  circular  in  the  absence  of  H will  remain,  after  H is  introduced, 
a precession  or  a nutation  (actually,  it  will  develop  into  an  epicycle). 

In  the  presence  of  E the  nutation  and  the  precession  are  no  longer 
circles,  but  spirals.  In  particular,  in  the  linear  theory  we  should  look 
for  those  spirals  which  are  “steady®;  naturally,  this  means  unchanging 
ratio  0/5.  This  can  indeed  be  the  case  if  the  torques  shrink  Or  stretch 
the  quasi-^arecession  in  the  same  manner,  proportionally,  as  the  angular 
velocity  T)  shrinks  or  stretches  the  yaw  6.  Two  such  possible  cases  are 
shown  in  Figure  ?•  A shrinking  spiral  (SE  regime)  means  1 lagging  behind 
the  plane  of  yaw;  then  M has  a component  which  shrinks  the  quasi-precession 
working  against  the  stretching  action  of  H.  An  expanding  spiral  (EE  regime 
means  L ahead  of  the  plane  of  yaw,  with  a component  of  M stretching  the 
quasi-nutation  (against  the  shrinking  action  of  H)  as  well  as  the  quasi- 
precession. 


* 


However,  H has  nothing  to  do  with  the  viscosity  of  the  air 
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Specifically,  our  requirement  now  is  that  the  proportions  of  the 
triangle  T^LA^  (rather  than  orOy  the  ratio  0/6)  remains  constant $ i.e., 

that  the  angular  rates  of  both  L and  A^  about  be  the  same,  and  the 

angle  between  quasi-precession  and  the  total  torque  be  the  same  as  the 
angle  between  the  plane  of  yaw  and  the  arrow  of  ,T)  . The  spiral  then, 
in  particular,  is  a logarithmic  spiral  (this  corresponds  exactly  to  the 
situation  in  the  mathematical  approach,  where  we  solve  a differential 
equation  by  inquiring  whether  there  may  not  be  solutions  the  Torn 
Xt 

e ).  In  particular,  we  may  assume  that  the  angles  UT.^  ^1^1  are 

small  (the  assumption  of  small  torques,  and  hence  small  damping  rates, 
is  - together  with  the  assumption  of  small  yaw  - a characteristic  feature) 
of  linear  theory).  Then  on  Figure  9 we  must  have 

(n/2  - p)  + U/&  = n/2  - y (8) 

Substituting  H *»  hr),  y)*  (A6n/B)0,  M * p6,  aid  f5  ^ -y0/(6-0)  we  may 
“ solve  for  y in  terms  of  the  ratio^/Bs 

y * h(Aoi^/Bn) (0/5) (l  - 0/6)/(20/5  - l)  (80 

Since  the  triangles  are  very  flat  (spiral  is  much  like  a circle), 
we  may  take  it  that  0/6  is  as  given  by  (3B);  then 

y » h^A^Cl  - ff2)/a  (8") 

From  the  definition  of  o in  (3B)  it  may  be  readily  noted  that 

1 - O2  m 

so  that  upon  substituting  this  into  (8n)  and  cancelling,  we  have  the 
simple  and  curious  result, 

y * ♦ h/Aco^c  (8A) 

i.e.,  the  spirals  of  nutation  and  precession  are  mirror  images  of  each 
other.  However,  since  nutation  goes  around  faster,  it  shrinks  faster. 
Multiplying  (BA)  by  the  angular  rates  as  given  by  (JO,  we  have  the  propor- 
tional radial  rates  (the  radial  components  of  T\  divided  by  6),  with  respect 
to  time,  as 

rrj/6  - + (h/AcOjff) (Aoo^/2B) (l+c)  * (h/2B)  (l*  l/ff)  (8B) 

To  compare  this  with  the  result  of  the  standard  "ballistic  texts  (cf 

li  2 2 

particularly,  EEL  1*1*6),  we  may  substitute  h » Kgpd4u  and  B - mk  d j then 
the  factor  h/B  is  (K^pcP /m)u/k2 d,  where  the  group  of  factors  in  the 


parentheses  is  usually  written  as  Jg,  so  that  (8B)  becomes  ( Jg/2k2 ) (u/d)  (l±l/o)  * 

The  factor  u/d  drops  out  when  we  change  the  independent  variable  from  time  to 
the  distance  travelled  by  the  shell  and  expressed  in  calibers.  Thus,  finally, 
the  proportional  radial  rates  are 

(Jj/2k2)(lil/a),  (BC) 

in  agreement  with  the  linear  theory. 

A curious  student  might  inquire  here  whether  some  of  our  simplifying 
assumptions  might  not  be  rescinded  at  this  point.  Indeed,  the  recission 
of  the  assumption  that  the  damping  is  small  would  be  a straight-forward  - 
if  laborious  - algebraic  task;  but  the  cumbersome  result  of  this  task  would 
not  be  particularly  useful,  simply  because  this  assumption  is,  generally, 
in  excellent  agreement  with  experiment.  More  important  - and  more  difficult  - 
would  be  a recission.  of  the  linearity.  On  a sphere  there  seems  to  be  no 
clean-cut  analog  of  the  logarithmic  spiral'9'  - or,  more  exactly,  it  is  hard 
to  choose  between  many  imperfect  analogs.  A crude  sketch  of  several  of  such 
analogs  is  given  in  Fig,  10.  We  shall  not  inquire  which  ones  of  these  analogs 
constitute  "possible"  yawing  motions,  and  under  which  assumptions;  it  seems 
probable,  for  instance,  that  several  of  such-like  motions  might  result  under 
the  identical  assumptions  about  M(6)  and  H(7}),  if  the  initial  conditions 
are  suitably  varied  - in  some  vague  analogy  to  the  blending  between  the 
nutation  and  the  precession  in  the  non-linear  regime.  What  matters  for  our 
present  purposes  is  that  the  idea  of  nutation  and  precession  - which  in  the 
non-linear  regime  has  already  lost  its  usefulness  of  superposition,  and 
survived  only  for  undamped  circular  motions  - becomes  even  more  complicated 
(viz.,  less  useful)  when  damping  is  introduced.  It  is  hoped  that  mathema- 
ticians will  eventually  clarify  these  concepts. 

So  far  we  have  only  retraced,  verbosely  and  laboriously,  some  of  the 
ground  wnich  is  explored  more  swiftly  and  elegantly  by  mathematics;  and  we 
have  refrained  from  doing  that  which  we  have,  in  effect,  promised  the  reader: 
to  defend  the  use,  in  ballistics,  of  that  obvious  mental  process  which  - for 
want  of  a better  name  - we  called  "qualitative  reasoning".  Let  us  attempt 
it  now.  As  soon  as  we  broaden  our  approach  by  concerning  ourselves  with 
an  approximate  (as  well  as  exact)  prediction,  the  non-linearity  ceases  to  be 
the  bugaboo  it  was  in  any  quantitative  reasoning. 

Let  us  consider  a typical  looping  epicycle  (linear  or  non-linear),  with 
loops  falling  short  of  the  origin.  We  sense  that  this  is  a case  in  *faich 
the  quasi -precession  predominates,  somehow,  over  the  quasi-nutation.  The 
phase  relation  between  A^,  L and  alternates  rapidly,  and  we  surmise  that 

the  alternating  stretching  and  shrinking  of  quasi-precession  and  quasi -nutation 
by  M,  on  the  whole,  cancel;  we  may  consider  then  only  the  direct  effects  of  H. 


Since  the  sum  of  angles  of  a spherical  triangle  depends  upon  its  size, 
there  can  be,  generally,  no  similar  triangles  (as  LT^)  on  the  sphere. 


bk 
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We  infer  that  the  quasi -nutation  definitely  damps.  As  to  the  quasi- 
precession, we  observe  that  it  is  stretched  when  the  angular  velocity  V\ 
is  large  (near  the  maximum  of  yaw) , and  shrunk  when  r\  is  small  (near  the 
minimum  of  yaw) ; we  infer  that  quasi-precession  grows',  but  only  slowly* 

We  therefore  infer  that  the  long-range  stability  is  present,  with  the 
following  qualification.  As  the  loops  shrink,  the  shrinkage  of  quasi- 
nutation slows  down,  and  the  growth  of  quasi-precession  becomes  more  steady; 
eventually,  we  no  longer  can  assume  that  the  phase  relations  of  A^  and  L 

cancel  out,  but  surmise  that  some  steady  phase  relation,  such  as  that  of 
of  Fig.  9b,  begins  to  be  felt,  i.e.,  even  the  quasi-nutation  will  begin  to 
grow;  i.e.,  the  long-range  stability  will  eventually  cease;  although,  by 
(8C),  the  growth  of  precession  will  be  much  slower  than  the  shrinkage  of 
the  nutation  had  been  (this  circumstance  is  also  obvious  from  a comparison 
of  Fig.  9a  and  9b:  in  precession  H is  smaller,  and  also  it  opposes,  rather 
than  adds,  in  some  ways,  the  action  of  M).  Now,  this  is  the  case  even  when 
we  start  with  a motion  in  which  the  quasi-precession  seems  to  predominate 
from  the  start.  Had  we  started  with  an  epicycle  looping  over  the  origin, 
the  long-range  stability  would  have  been  even  more  pronounced,  since  the 
situation  would  more  often  resemble  Fig.  9a;  yet  we  cannot  claim  the  long- 
range  stability  for  an  indefinite  range,  as  long  as  the  situation  is  sot 
precisely  that  of  Fig.  9a  - for  a perturbation  of  the  steady  phase  relations 
of  Fig*  9a  would  tend  to  cause  a growth  of  quasi -process  ion  precisely  at 
those  times  when  7}  is  greater  * 

Thus,  our  qualitative  conclusion  (without  going  into  any  mathematics ) 
is  that  quasi-nutation  damps  strongly  when  it  is  a sizable  component  of 
yaw,  grows  slowly  only  when  it  is  an  unimportant  component;  the  quasi - 
precession  damps  strongly  when  it  is  an  unimportant  component  of  yaw,  grows 
slowly  when  it  is  the  principal  component  of  yaw*  This  conclusion  falls 
pretty  close  to  what  is  predicted  by  the  exact  linear  theory;  yet,  it  is 
completely  free  from  the  assumptions  of  linearity. 

MAGNUS  TORQUE 

A positive  Magnus  torque*  on  all  our  sketches  pulls  L up.  Thus  in  a 
general  epicycle  (if  the  oscillatory  effects  of  M may  be  presumed  to  cancel 
out  on  the  whole)  T builds  up  the  quasi -precession.  The  quasi-nutation  is 
on  the  whole  shrunk:  for  it  is  shrunk  when  5 and  P are  large,  stretched 
only  when  6 and  T are  small.  The  steady  spiral  motions . (Fig.  11)  are  possi- 
ble if 


(ji/2  - fl)  + T/fc  - n/2  - T ( 9 ) 

Substituting  p and  M as  before,  and  T « tS,  we  can  solve  for  yt 
T - (t/u)(l  - 0/6)/(2j2f/6  - 1)  (90 


4 

Magnus  force  is  the  force  which  deflects  a spinning  ball,  and  Magnus 
torque  is  the  torque  produced  by  the  Magnus  force*  The  reader  is 
expected  to  be  familiar  with  this  (and  other)  historical  background 
of  Magnus  effects. 
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Using  (3A)  again,  we  have 

Y • t (t/2n)(l  i a) /a  • - (t/2n)  (l  ? 1/ct),  (9*) 


i.e.,  the  two  spirals  are  no  longer  mirror  images » the  precession  spiral 
is  steeper*  Finally,  substituting  p.  ■ A2aa^2(l  - j2)/UB,  t * and 

A ■ and  reverting  to  the  notation  of  BRL  U*6,  we  have  the  proportional 

Si 

rates  of  growth  of  nutation  and  precession  as 

♦ (yia*)M 

in  agreement  with  the  linear  theory. 


(9C) 


The  effect  of  T is  somewhat  similar  to  that  of  H,  the  differences  lying 
in  the  manner  of  dependence  of  damping  upon  the  stability  parameter  a,  and 
also  in  the  fact  that  T can  be  negative  as  well  as  positive*  Many  interesting 
examples  of  the  interplay  of  these  facts  are  discussed  extensively  elsewhere 
(cf.,  e.g.,  BRL  668  and  BRLM  682).  It  should  be  noted  that  the  study  of  the 
effects  of  T seems  to  be  pne  occasion  where  it  is  desirable  - because  of  the 
recently-discovered  strong  non-linearity  of  T - to  substitute  for  the  lineari- 
zation (9')  an  exact,  non-linear,  solution  of  (9). 

An  illustration  of  the  effects  of  H and  X (as  well  as  the  effect  of  the 
aerodynamic  forces,  which  we  shall  discuss  presently)  is  provided  by  the 
spinning  top*  The  friction  at  the  tip  of  the  top  acts  as  some  analog  of 
aerodynamic  cross-forces,  and  its  torque  as  some  analog  of  certain  aerodynamic 
torques*  The  analogy  is  imperfect  because  the  laws  of  Coulomb  friction  are 
not  the  same  as  those  of  the  aerodynamic  forces  and  torques  (the  friction 
force  depends  only  upon  the  direction,  rather  than  the  magnitude,  of  the 
velocity  at  which  the  tip  of  the  top  rubs  against  the  table).  The  greater 
complexity  of  the  top's  motion  (in  comparison  with  that  of  the  shell)  results 
in  the  apparent  lively  and  capricious  behaviour  of  the  top;  but  there  can  be 
no  doubt  that  some  of  the  dispersion  of  artillery  fire  is  akin  to  the  wandering 
of  the  top. 


At  the  start  the  center  of  gravity  of  the  top  is  substantially  stationary, 
while  the  tip  rapidly  swings  out,  sliding,  substantially  in  a circle*  The 
frictional  force  is  substantially  perpendicular  to  the  plane  of  yaw,  and  acts 
as  a positive  Synge  force  S (or  as  a positive  Magnus  force  F)j  it  produces 
some  analogs  of  H (or  a positive  T)  - and  also  a slight  spin-accelerating 
torque*  As  a result,  the  nutation  damps  out  (rather  abruptly),  and  presently 
only  the  pure  precession  remains*  This  precession  may  (because  of  H and 
positive  T)  grow  for  a while j but  the  damping  is  presently  taken  over  by  other 
processes*  Being  (in  the  slow  precession)  rather  steady  in  its  direction,  the 
frictional  force  now  causes  same  horizontal  velocity  of  the  c*g.  across  the 
plane  of  yaw:  this  tends  to  neutralize  the  rubbing  velocity  perpendicular 
to  the  plane  of  yaw,  so  that  the  tip  tends  to  start  rolling,  rather  than  sliding, 
on  the  surface  of  the  table.  However,  as  the  plane  of  yaw  continues  processing, 
the  rubbing  velocity  develops  a component  in  the  plane  of  yaw  (the  top  starts, 
literally,  to  "draw*  on  the  surface).  The  frictional  force  thus  develops  a 
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component  analogous  to  the  lift  L,  which  (as  it  will  be  presently  shown)  tends 
to  damp  the  precession  (there  being  no  longer  a nutation  to  grow) . . It  seems 
to  be  an  experimental  fact  (undoubtedly  provable  by  analysis)  that  the  effect 
o t L overcomes  that  of  H and  T.  This  "dragging11  neutralizes  the  velocity  of 
the  c.g.,  and  finally  the  vector  a comes  within  the  radius  of  the  tip;  H (or  T) 
become  negative,  and  the  top  rather  suddenly  "goes  to  sleep". 

Particularly  interesting  is  the  kind  of  the  top  which  turns  upside  down. 

Such  a top  starts  as  statically  stables  viz.,  the  radius  of  its  tip  is  large, 
and  the  metacenter*  of  the  top  is  above  the  center  of  gravity  (rather  than  at 
the  tip).  For  the  same  reason  the  rubbing  velocity  due  to  spin  is  generally 
greater  than  the  rubbing  velocity  due  to  yawing;  the  friction  is  analogous, 
therefore,  to  a negative  Magnus  torque.  Now  it  is  the  precession  which  damps 
out  rapidly,  while  the  nutation  builds  up  to  n/2  and  beyond.  While  the  vector 
of  angular  velocity  remains  substantially  vertical,  it  turns  with  respect  to 
the  top;  for  an  instant  the  top  is  in  the  pure  cartwheeling  motion,  without 
any  axial  spin.  Thereafter  the  motion  can  be  considered  as  the  nutation  of 
a statically— uns table  top,  with  a positive  Magnus  torque  - which  presently 
commences  damping.  This  is  as  though  a fin-stabilized  shell  turned  around 
and  flew  fins  forward;  in  fact,  the  early  stages  of  this  motion  are  not  unlike 
the  motion  of  a mortar  shell  landing  "short". 

OTHER  CROSS-TORQUES 

The  cross-torques  M,  H and  T pretty  well  complete  the  cross-torque  part 
of  Kelley-McShane  matrix  (7);  symmetry  demands  only  that  we  mention  the  Magnus - 
cross-torque  XT.  This  is  one  of  the  marginally  relevant  components,  given 
rather  an  undue  prominence  by  the  Maple-Synge  symmetry  arguments;  it  now  seems, 
for  instance,  that  various  non-linear  components'"* of  T are  more  important  than 
XT  (it  should  be  noted,  however,  that  XT  has  not  been  measured  at  really  large 
yaws).  On  our  sketches  XT  pulls  L across  the  quasi-nutation;  thus,  in  the 
general  epicycle  its  damping  effects  rather  cancel,  and  in  circular  yawing 
motions  it  is  rather  difficult  (through  not  impossible***)  to  distinguish  its 
effect  from  that  of  M,  since  it  affects  only  the  angular  rates.  This  is  in 
agreement  with  the  theory  of  BRL  UU6,  which  shows  that  XT  produces  only  a 
negligible  effect  on  the  angular  rates. 

We  might  remark  on  the  apparent  (through  so  far  not  established)  checker- 
board and  hierarchical  character  of  this  matrix.  Die  linear  M can  be  generalized 
into  a function  of  both  yaw  and  spin,  odd  in  yaw  and  even  in  spin;  all  terms 
of  such  an  expansion  of  M seem  to  affect  only  the  angular  rates,  and  leave 
damping  unaffected.  The  torque  T,  odd  in  both  yaw  .and  spin  (and  perpendicular, 
by  Maple-Synge  theory,  to  M) , appears  as  a sort  of  refinement  on  M;  all  its 

• Which  is  a point  on  the  axis  of  the  top  (or  fixed  with  respect  to  a ship) 
thru  which  the  reaction  of  the  table  (or  of  water)  acts,  for  small  angles 
of  yaw. 

#*  Cf  BRLM-682,  "Spiral  Tawing  Motions  of  31  mm  Shell  M£6:  a Study  in  Non- 
Linear  Theory"  (1953). 

-**#  Cf  BRL  882.  - 
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terns  seem  not  to  affect  angular  rates,  but  are  very  essential  iir  damping. 

How,  the  angular  velocity Y)  in  itself  is  a sort  of  refinement  on  yaw  6;  thus, 

H can  be  accounted  for*  by  the  variation  of  the  instantaneous  angle  of  attack 
along  the  axis  of  a yawing  shell.  Thus  R,  analogously  to  T,  seems  a refinement 
on  M;  it  can  be  an  even  function  of  spin,  does  not  affect  the  angular  rates, 
and  affects  damping.  Then  XT  appears  a refinement  on  both  B and  Ts  it  stands 
in  the  same  relation  to  H as  T to  M,  and  to  T as  H to  H;  it  is  odd  both  in  4) 
and  spin,  affects  angular  rates  and  does  not  affect  damping.  This  is  as  though 
the  linear  theory  elevated  XT  on  the  principle  of  primogeniture,  and  disinherited 
its  more  important  cousins,  e.g.,  non-linearities  of  T. 

CROSS-FORCES 

The  cross-forces  - vis.,  the  aerodynamic  forces  N,  5,  F and  XF  of  matrix 
(7),  as  well  as  the  component  of  gravity  normal  to  the  trajectory  - affect  the 
yawing  motion  by  producing  a swerve,  l.e.,  by  moving  point  T^,  or  by  shifting 

the  base  of  the  "vector*  of  quasi-precession,  rather  than  that  of  the  quasi- 
nutation.  This  is  considerably  different  from  the  effect  of  the  cross-torques , 
and  our  "qualitative  reasoning",  accordingly,  will  have  to  beTconsiderably 
modified.  Also,  it  will  not  be  satisfying  as  it  was  with  the  cros s -torques • 

The  difficulty  is  quite  real:  we  are  dealing  with  a more  complicated  phenomenon 
than  a mere  angular  motion,  viz.,  with  the  effect  of  swerve  upon  the  yawing 
motion.  Mathematically,  this  is  the  step  of  coupling  together  two  separate 
differential  equations  - a step  laborious  even  in  the  linear .theory.  Indicative 
of  this  distinction  is  the  obvious  "cumulative*  character  of  the  effect  of 
motion  of  T^  upon  the  motion  of  L:  the  velocity  of  must  have  existed  for 

some  length  of  time,  and  in  a substantially  the  same  direction,  before  its 
effect  is  felt;  or,  the  position  of  T^  in  the  configuration  of  the  triangle 

T^  LA1  reflects  not  the  instantaneous  velocity  of  T^,  but  that  velocity  which 

had  existed  same  time  previously. 

Lift.  In  advanced  exterior-ballistics  texts  (e*g»,  ERL  Ui6)  the  forces 
in  the  plane  of  yaw  are  resolved  usually  into  the  normal  force  K (perpendicular 
to  the  axis  of  the  shell)  and  axial  drag  DA.  It  seams  more  natural  to  resolve 
them  into  lift  £ (perpendicular  to  the  trajectory)  and  drag  D (opposed  to  T^); 

in  our  case  this  is  also  more  convenient. 

On  all  our  sketches  lift  pulls  T^  up.  To  the  first  approximation  - if  we 

can  disregard  the  effect  of  M on  the  damping  of  the  yawing' motion  - it  is 
obvious  that  quasi-precession  is  shrunk  by  the  lift  (those  rare  occasions  when 
L on  our  sketches  is  below  T^  represent  minima  oi,  yaw,  when  lift  is  small). 

The  effect  of  lift  on  the  quasi -nutation  is  less  obvious,  and  in  fact  is  less 
important.  To  the  first  approximation  we  could  simply  say  that  a motion  of  T^ 

has  no  effect  upon  L and  A, , i.e.,  has  no  effect  on  the  quasi-nutation.  Indeed, 
that  would  not  be  very  far  from  the  truth. 


Cf  BRL  668;  also,  recent  HRL  monograph,  "Exact  Re-Statement  of  the  Equations 
of  Motion  of  a Shell". 
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On  a closer  inspection  we  observe  that  the  motion  of  does  affect  L 

by  changing  the  yaw,  and  therefore  changing  the  torque  M*  Specifically,  we 
may  distinguish  two  effects  of  the  motion  of  1^*  ” 

The  first  one  is  as  follows.  In  its  continuous  running  toward  A^,  the 

point  T^  fails,  so  to  say,  to  "lead11  A^.  The  plane  of  yaw,  therefore,  turns 

more  rapidly,  and  so  does  the  arrow  of  M on  our  sketches.  Inspecting  our 
sketches  (either  one,  or  every  one,  according  to  the  reader's  patience)  we 
observe  that  the  faster  turning  of  M tends  to  pull  L away  from  A^;  we  conclude 

that  lift  tends  to  stretch  the  quasi-nutation. 

We  are  stepping  now  close  to  a pitfall.  For  instance,  the  linear  theory 
shows  that  precession  and  nutation  (under  the  assumptions  which  we  have  so 
far  accepted,  viz.,  a constant  velocity  and  spin  of  the  shell). have  constant 
proportions  0/5;  but  precession  is  damped  by  lift,  while  nutation  is  built 
up.  Hence  the  quasi-precession  in  nutation  is  stretched,  and  quasi -nutation 
in  precession  is  shrunk  - in  apparent  contradiction  to  our  preliminary  obser- 
vations . The  answer  lies,  again,  in  the  neglect  of  the  effect  of  M,  and  in 
the  unwarranted  assumption  that  a motion  representing  a certain  proportion 
of  precession  and  nutation  without  the  lift  will  retain  this  proportion  after 
the  lift  is  introduced.  It  is  easy  to  confuse  two  different  problems  : one 
(ours),  how  does  the  lift  affect  a motion  with  given  initial  conditions;  and 
the  other  (practically  habitual  in  the  milieu  of  mathematical  approach),  how 
does  a possible  spiral  motion  in  the  presence  of  lift  differ  from  a possible 
spiral  motion  without  lift. 

We  have  an  inkling  of  this  in  considering  the  second  effect:  as  M turns, 
it  also  (on  the  whole)  tends  to  diminish  in  magnitude,  because  of  the  gener- 
ally diminishing  yatw.  The  angular  rate  of  quasi-precession,  we  might  judge, 
is  not  affected  much,  for  the  magnitude  of  the  quasi-precession  also  diminishes. 
Bnt,  since  the  velocity  of  is  independent  of  the  position  of  T^,  the  angular 

rats  of  about  T^  is,  on  the  whole,  increased;  i.e.,  the  synchronism  of  L 

and  A^a  dr  the  spiral  character  of  the  motion,  will  not  he  generally  preserved  - 

unless  the  phase  relations  of  the  triangle  LA^  are  such  that  the  magnitude 

of  quasi-nutation  is  changing  in  some  particular  manner.  In  this  way  we  come 
to  consider  the  spiral  motions  again.  Two  possible  configurations  are  shown 
on  Fig*  13,  which  explains  the  physical  mechanism  of  the  well-known  result  of 
the  linear  theory.  Thus,  in  precession  (Fig.  13b)  the  quasi-precession  is 
shrunk  not  only  by  the  lift,  but  also  by  the  M-torque;  the  quasi -nutation  is 
simultaneously  strongly  shrunk  by  the  M-torque.  In  nutation  (Fig.  13a)  the 
M-torque  stretches  the  quasi-precession  more  strongly  that  it  is  shrunk  by 
the  lift;  and  simultaneously,  the  M-torque  slightly  stretches  the  quasi- 
nutation. 

Our  qualitative  conclusion  (completely  free  from  the  assumptions  of 
linearity)  is,  therefore,  as  follows.  The  quasi-precession  is  damped  by  lift 
strongly  if  it  i3  an  essential  component  of  yaw,  grows  slightly  if  it  is  an 
unimportant  component.  The  quasi-nutation  grows  slightly  if  it  is  the  princi- 
pal component  of  yaw,  shrinks  strongly  if  it  is  an  unimportant  component. 
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F,g  !3:  Effect  of  Lift 


Fig  14:  Yaw  of  Re. pose,  in  the  Presence 
of  Damping  Torque  (an  exainp/e.) 


The  effect  of  lift  should  be  juxtaposed  with  that  of  H-torque.  In  fact, 
these  two  refinements  on  an  undamped  epicyclic  motion  form  a would-be  satis- 
factory system  for  damping  of  yawing  motion  in  the  absence  of  Magnus  effects: 

One  damps  nutation,  the  other  - precession.  This  fact  has  been  brought  out 
particularly  clearly  by  Davis  and  Foilin'*,  to  whose  work  I take  pleasure  to 
acknowledge  my  indebtedness  once  again.  The  student  of  ballistics  should  be 
reminded  at  this  point  that  the  Magnus  effects  had  to  be  introduced,  historical- 
ly, precisely  because  that  system  failed  to  account  satisfactorily  for  the 
totality  of  experimental  information. 

Let  us  now  revert  to  the  quantitative  evaluation  of  the  effect  of  lift 
in  the  linear  case.  Instead  of  considering  the  angles,  as  we  did  in  (8)  and 
(9),  we  should  now  consider  (cf*  Fig.  13)  the  uroportional  radial  rates,  i.e., 
idle  rates  of  growth  divided  by  the  magnitudes.  To  preserve  similarity, 

(rr)- -C/*u)/&  - (pMA  /mu)/(6  - 0)  (10) 

Substituting  Z - /6,  M ■ |a6,  L ■ Aos^,  T]  » (AoOj/B)0,  and  p » y0/(5  - 0) 
and  solving  for  y,  we  have 

T - - (jf/mu)  ( B/Aos^)  (l  - 0/$)/  [(1  - 0/5)2  - WA2«£  ] (10*) 

Substituting,  from  (3*)  and  (3B),  Bn/A2ea^  * (1  - a2) A and  requiring  that 

the  ratio  0/&  remain  constant  as  in  (3A) , we  have  the  simple  result,  somewhat 
analogous  to  (8A): 

t * i (//m)  (B/Aa^)/a  (10A) 

In  the  presence  of  swerve  we  no  longer  can  compute  the  proportional  radial 
rates  (with  respect  to  time)  as  t\y/&  - which  we  did  in  (8B) } nor  can  we  say  - 
in  analogy  to  (8A)  - that  the  two  spirals  are  mirror  images  of  each  other; 
rather,  these  rates  are  given  by  the  left-hand  side  of  (10).  Substituting 
into  (10)  / ■ J^pd2^,  and  reverting  to  the  notation  of  BEL  1*1*6,  we  have  the 

proportional  radial  rates  with  respect  to  travel  as 

- (Jl/2)(1  + 1/0),  (IOC) 

in  agreement  with  the  linear  theory. 

There  is  some  resemblance  between  (IOC)  and  (8C),  but  it  should  not 
distract  our  attention  from  the  fact  that  the  mechanisms  of  the  spirals  of 
(IOC)  and  (8C)  are  essentially  different.  The  spirals  of  (8C)  were  yawing 
motions  about  a straight-line  trajectory,  those  of  (lOC)  represent  an  effect 
of  swerve.  Similarly,  some  further  refinements  on  (l£))  are  possible,  but  do 
not  seem  especially  profitable. 


* OSRD  Report  2 #9,  CIT  JPC  29. 
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Magnus  Force,  While  the  appreciation  of  the  Magnus  effects  arose 
originally  from  the  knowledge  of  the  -Existence  of  the  Magnus  force  F (cf . 
the  swerve  of  a spinning  ball,  Magnus's  wind-tunnel  experiments  with  spinning 
cylinders  in  cross-flow,  Flettner's  rotor,  etc.),  with  a shell  the  most  im- 
portant effect  of  Magnus  force  F is  that  it  produces  the  Magnus  torque  T - 
whose  effects  have  already  been  discussed.  Magnus  force  proper  does  produce 
some  legitimate  swerve*,  (and  in  a spinning  top,  the  friction  at  the  tip 
produces  some  wandering  of  the  top) ; but  its  direct  effect  on  the  yawing  motion 
per  se  is,  fortunately,  negligible.  A positive  F on  all  our  sketches  pulls 
to  the  left?  therefore,  it  affects  mostly  the  angular  rates  (slightly 

increasing  them),  rather  than  the  damping  rates  - and  its  effect  is  difficult 
to  distinguish  from  that  of  a change  in  H**. 

Synge  Force.  The  force  S which  produces  the  damping  torque  H has  been 
introduce!  into  ballistics  for  the  sake  of  a certain  formal  mathematical 
consistency.  This  force  is  closely  analogous  to  Magnus  farce.  On  all  our 
sketches  this  force  pulls  across  the  direction  of  the  "vector*  of  quasi- 
nutation I Thus  in  the  general  epi cyclic  motion  its  effects  rather  cancel  out 
(even  more  so  than  was  the  case  with  F),  and  in  the  spiral  motions  it  affects 
only  the  angular  rates  (rather  than  the  damping  rates);  It  does  have  some 
effect  on  the  swerve,  but  in  the  study  of  small  yawing  motion  it  can  be  dis- 
regarded. 

Our  qualitative  conclusions  about  the  unimportance  of  F and  S are  in 
agreement  with  the  results  of  the  linear  theory. 

Cross-Magnus  Force.  The  force  XF  ifaich  produces  the  Magnus  cross-torque 
XT  needs  mentioning  only  for  the  sake  of  completing  the  linear  matrix  (7);  it 
survives  in  the  linearisation  of  the  Maple-Synge  theory  only  for  the  reason 
that  Synge's  aerodynamic  hypothesis  treats  the  yawing  velocity T)  as  being  on 
par  with  the  yaw  6 (rather  than  as  a refinement  on  6,  as  we  hold).  In  the 
results  of  the  linear  theory  of  yawing  motion  the  effect  of  this  force  is 
apparently  negligible.  Yet  we  may  surmise  that  the  cross -force  part  of  the 
matrix  (7)  may  possess  the  same  checker-board  character  which  we  have  detected 
in  the  cross-torque  part  of  that  matrix;  then  IF  should  affect  damping.  On 
our  sketches  XF  pulls  along  the  direction  of  quasi -nutation;  90  that,  indeed, 

it  affects  damping  in  the  basic,  spiral,  motions.  It  is  a simple  matter  to 
investigate  the  effect  of  XF  on  spiral  motions  qualitatively,  with  the  assump- 
tions of  linearity  - analogously  to  (10) . The  result  is  that  the  radial  pro- 
portional rates  of  growth  with  respect  to  the  distance  in  calibers  are 

- JnVkaa  Cue) 


* Cf.,  e.g.,  ERL  703,  "On  Jump  Due  to  Mu2zle  Disturbances";  also  BRLM  682. 

An  evaluation  of  the  effect  of  F on  the  angular  rates  could  readily  be  made, 
but  does  not  appear  profitable?  for,  with  the  assumptions  usually  made  in 
the  linear  theory,  it  turns  out  that  this  effect  is  negligible. 
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This  is  in  agreement  with  what  the  result  of  linear  theory  would,  have  teen 
if  terns  of  second  order  in  JTs  had  been  kept.  Thus,  XF  does  not  matter  in 
yawing  motion  simply  because  its  effects  are  small.  We  may  note  that  our 
surmise  about  the  checker-board  character  of  (7)  is  confirmed,  but  is  over- 
whelmed by  the  hierarchical  character  of  that  matrix*  If  the  reader  will 
tolerate  another  figure  of  speech,  the  primogeniture  claims  of  XF,  filed  with 
the  Maple -Synge  theory,  are  rejected  even  by  the  Kelley-McShane  theory  which 
is  descended  from  Nielsen-Maple -Synge  concepts. 

The  checker-board  character  of  (7)  cun  be  summed  up  by  writing 


M(a) 

T(d) 

N(d) 

and 

F(a) 

H(d) 

XT  (a) 

S(a) 

XF(d) , 

where  (a)  and  (d)  stands  for  "affects  angular  rates*  and  “affects  damping 
rates".  The  hierarchical  character  of  (7)  can  be  summed  up  by  ranking  its 
components  in  the  order  of  their  importance: 


1st  rank 
2nd  * 
3rd  * 
Uth  * 


M 

(a) 

determines  the  basic  character  of  motion 

H,  T, 

N 

(d) 

determine  damping 

XT,  F, 

S 

(a) 

refinement  on  angular  rates 

XF 

(d) 

refinement  on  damping  rates. 

The  usefulness  of  this  trivial  classification  lies  in  the  fact  that  it 
shows  up  certain  weaknesses  of  flielsen-Synge-Kelley-McShane  linearization  of 
Maple-Synge  theory,  and  helps  to  organize  the  methods  of  further  attack:  thus, 
if  we  study  the  damping  rates,  no  direct  effect  of  XT,  F and  S can  be  expected, 
while  refinements  on  H,  T and  N are  just  as  legitimate  - and  may  be  more 
important  - than  the  measurement  of  the  "first-order*  term  XF. 


TAW  OF  REPOSE 

The  effect  of  the  force  of  gravity  acting  across  the  trajectory,  which 
causes  the  curvature  (as  distinguished  from  the  swerve)  of  the  trajectory,  is 
considerably  different  from  that  of  the  aerodynamic  cross-forces  - for  the 
simple  reason  that  this  force  is  relatively  constant  and  independent  of  yaw 
and  yawing  velocity.  Mathematically  this  force  manifests'  itself  an  inhomo- 
geneity of  the  differential  equation  (rather  by  a change  in  the  coefficients 
of  the  equation)  and  results  in  the  fact  that  the  general  (linear)  solution 
consisting  of  precession  and  nutation  must  be  augmented  by  the  third  component, 
called  yaw  of  repose. 


We  may  arrive  at  the  same  conclusion  by  an  argument  analogous  to  that  by 
tfiich  we  have  introduced  the  generalized  (non-linear)  precession  and  nutation. 
The  point  "falls",  in  the  surface  of  our  sphere,  (pf.  Fig.  Hi)  at  the 

velocity  gcosS/u  (where  g is  the  acceleration  of  the  gravity  and  & la  the  angle 
of  inclination  of  the  trajectory).  If  the  shell  is  to  trail  properly,  the 
vectorial  velocity  of  on  thi3  sphere  must  be  approximately  the  same.  Drop- 
ping the  convention  of  drawing  above  T^,  we  observe  that  L must  be  directly 
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to  the  left  of  i.e.,  there  must  exist  a horizontal  quasi-nutation,  pointing 

(from  L to  A^)  to  the  rights  Moreover,  the  point  L of  this  sphere  must  have 

substantially  the  same  vectorial  velocity,  i.e.,  there  must  exist  suitable 
torques j vti  assume  these  to  consist  largely  of  the  H- torque . Then  there  must 
exist  a yaw  (from  to  A^  pointing  directly  to  the  right.  Thus  T^,  L and 

are  co -planar  and  on  our  sketches  the  situation  appears  as  a circular  yaw, 

except  that  the  motion  of  causes  the  triplet  of  these  points  to  move  down 

without  any  rotation  noticeable  in  our  sketches. 


The  condition  of  the  "steadiness*  of  this  motion  is 


gcosS/u  ■ M(6)/L  -/"lain# 

(12) 

which,  on  assumptions  of  linearity,  becomes  simply 

gcos6/u  » ^6/Ao^r  )0i 

(12 ») 

and  from  which  there  follow  the  simple  results 

0/6  - lAs 

(12") 

and 

6 - (Ao^/p.)  (gcosO/u) 

(12A) 

The  expression  (12A)  for  the  yaw  of  repose  may  be  rewritten  in  a number 
of  ways.  For  instance,  in  the  notation  similar  to  that  of  HRT-  UU6  and  BEL  668 
we  would  introduce  the  dimensionless  quantities  v “ oijd/u  » spin  per  caliber 

of  travel,  and  Jqc  ■ gdcosfi/u^ ; then  (l2A)  becomes 

■vV'V 

(12C) 

in  agreement  with  the  skeletonized  form  for  the  yaw  of  repose  in  the  linear 
theory. 

The  yaw  of  repose  is  considered,  naturally,  a basic  form  of  yawing  motion, 
analogous  to  the  precession  and  the  nutation  of  (3B).  Its  superposabillty  with 
the  nutation  and  the  precession  is,  of  course,  distinctly  a feature  of  the 
linear  theory.  In  attempting  to  combine  the  yaw  of  repose  with  the  epicycle 
in  the  general  non-linear  theory  we  would  meet,  unfortunately,  the  difficulties 
analogous  to  those  we  met  in  attempting  to  pass  from  the  linear  to  non-linear 
epicycle.  A study  of  the  non-linearities  of  the  yaw  of  repose  existing  alone 
(analogous  to  our  review  of  the  theory  of  the  spinning  top)  does  not  appear 
profitable. 

It  is  a simple  matter,  however,  to  generalize  (12),  within  the  assumptions 
of  linear  theory,  by  perturbing  it  for  the  effect  of  other  forces  and  torques. 
At  this  stage  we  may  readily  make  the  perturbations  for  the  cross-torques  and 
cross-forces:  these  would  be  very  similar  to  those  of  (8),  (9),  (10)  and  (11), 


except  that  the  conditions  of  spirality  would  be  replaced  by  the  conditions 
of  a parallel  synchronous  motion,  e.g,,  (3B)  replaced  by  (12*).  Just  to 
indicate  the  simplicity  of  such  processes,  we  might  make  the  perturbation  of 
(12‘)  for  the  H-torque.  From  Fig,  1U  it  is  obvious  that  we  must  have 

|j,6y  - h7]  (13) 

from  which,  using  (12"), 

Y « h/Ao^  (13A) 

and  the  downward  component  of  the  yaw  of  repose  can  be  shown  to  be 

6t  - WV  <>*> 

in  agreement  with  the  linear  theory. 

Additional  perturbations  on  the  yaw  of  repose  (including  the  perturbations 
for  the  axial  forces  and  torques,  which  we  have  not  made  as  yet  for  precession 
and  nutation)  we  shall  leave  as  an  exercise  for  the  reader:  and  we  need  hardly 
mention  that  our  method  can  reproduce  all  the  results  of  the  linear  theory, 

We  might  only  remark  that  in  any  respectable  trajectory  the  yaw  of  repose  is 
extremely  small,  and  is,  therefore,  in  itself  only  a minor  refinement  on  the 
epicycle.  Experience  shows,  in  some  cases,  that  an  apparently  large  summital 
yaw  is  very  likely  due  to  a survival  of  the  epicyclic  yaw,  i.e.,  insufficient 
damping.  We  might  also  mention  that  V unlike  the  aerodynamic  Jrs,  is  pro- 
portional to  the  caliber  - so  that  the  yaw  of  repose  must  be  scaled  from  one 
caliber  to  another  by  other  rules  than  those  of  the  epicycle*, 

FORCES  ALOMS  THE  TRAJECTORY 

Let  us  now  consider  the  effects  of  drag  (D)  and  of  the  component  of 
gravity  (mgsiofl)  acting  in  the  direction  opposite  to  the  trajectory  vector 
o Their  effects  differ  considerably  from  those  of  the  cross-forces 

(just  as  the  *WKB*  perturbation  of  BRL  kU6  differs  from  the  inclusion  of  the 
complex  components  into  the  coefficients  of  a differential  equation).  They 
affect  the  yawing  motion  by  changing  the  velocity  of  the  shell,  and  hence 
changing  the  overturning  moment  M and  the  stability  factor  s»  In  precession 
and  nutation  - cf  * (3B)  and  (IvB)  - the  proportions  of  the  spherical  triangle 
T^LA^  are  no  longer  constant,  even  in  the  linear  theory. 

Let  us  start  these  considerations  first  in  the  most  general  (and  that 
necessarily  means  vague)  way,  particularizing  them  presently  by  resorting 
to  the  spiral  motions,  and  finally,  to  the  assumptions  of  linearity*  Let 
us  start  at  the  instant  of  minimum  yaw  (either  W or  E sketch) . We  note  that 
throughout  the  regimes  of  increasing  yaw  (NW,  N and  NE  sketches)  the  motion 
of  the  point  L,  caused  by  M,  contributes  to  the  increase  of  yaw  - specifically. 


Cf.  E8LM  No,  833,  "Trajectory  Models  in  Mortar  Fire*. 
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by  tending  to  stretch  both  the  quasi-precession  and  the  quasi -nutation.  Hence 
the  decrease  of  M,  caused  by  the  loss  of  velocity*  decreases  the  growth  of  yaw. 

In  the  regimes  of  the  decreasing  yaw  (5E,  S and  SW  sketches)  the  motion 
of  L tends  to  shrink  both  the  quasi-precession  and  the  quasi-nutation.  This 
fact*  taken  alone*  does  not  give  us  the  information  as  the  total  effect  on 
yaws  both  the  quasi-precession  and  the  quasi -nutation  are  decreased  by  the 
loss  of  velocity  of  the  diell;  but  the  minimum  yaw  is  the  difference  between 
the  two.  We  can  observe  then  that  in  these  regimes  the  motion  of  L impedes 
the  shrinkage  of  yaw.  Hence  the  decrease  of  the  speed  of  L - by  decreasing 
the  impeding  of  the  shrinkage  of  yaw  - increases  the  shrinkage. 

We  conclude  tentatively  that  the  loss  of  velocity  always  damps  the  yaw 
of  a spin-stabilized  projectile. 

We  may  note  that  this  damping  is  associated  with  an  increase  of  the 
stability  factor;  and  it  is  interesting  that  our  intuitive  sense  of  the  ex- 
istence of  a relation  between  the  stability  and  the  amplitude  (or  damping) 
of  the  yaw  is*  in  a way*  confirmed  J* ..  We  may  also  note  that  at  this  stage 
we  are  not  able  as  yet  to  make  a definite  statement  about  the  shrinkage  of 
the  quasi -precession  and  the  quasi-nutation;  this  corresponds  to  the  fact 
that  we  made  no  committment  about  the  shape  of  the  epicycle*  viz.*  about  the 
proportions  of  the  precession  and  nutation  (or  quasi -precession  and  quasi- 
nutation) in  this  epicycle.  It  may  be  seen*  for  instance*  from  (3B)  that 
different  behavior  of  the  quasi-precession  and  quasi-nutation  may  be  expected* 
depending  upon  whether  the  nutation  or  the  precession  predominates. 

We  arrive  at  the  consideration  of  spiral  motions  by  a process  somewhat 
similar  to  that  we  used  in  discussing  the  effects  of  lift.  The  angular  velocity 
-A  of  quasi -nutation  about  the  instantaneous  position  of  point  L remains  con- 
stant, but  the  linear  velocity  of  point  L in  our  sphere  decreases!  thus  a 
possible  synchronism  of  the  quasi -precession  with  the  plane  of  yaw  tends  to 
be  disrupted.  let  this  synchronism  may  be- preserved,  approximately*  if  the 
* magnitude  of  the  quasi-precession  decreases*  t6o.  This  can  be  done  by  the 
torque  M if  L lags  the  plane  of  yaw;  the  situation  must  be  then  described  by 
a SE  sketch*  and  the  quasi-nutation  must  shrink*  too.  But  the  spiral  motions 
must  now  be  considered  in  conjunction  with  the  fact  that  the  ratio  J0/6  from 
(35)  or  (liB)  is  no  longer  a constant;  i.e»*  the  quasi-nutation  and  quasi-, 
precession  are  decreasing  at  different  proportional  rates*'  the  angles  p and 
Y in  the  triangle  T^LA^  are  changing*  and  the  synchronism  of ; the  quasi-  

precession  with  the  plane  of  yaw  can  no  longer  be  exact.  Specifically*  in 
nutation  (Fig.  l$a)  the  quasi -precession  will  shrink  more  rapidly*  and  the 


With  fin-stabilized  projectiles  the  loss  of  velocity  builds  up  the  yaw  - 
but  there  the  loss  of  velocity  mates  the  shell  more  sluggish  in  its  yawing 
motion,  or  - in  our  intuitive  sense  - less  stable. 
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ratio  0/6  will  increase,  as  it  should  by  (3B);  the  angle  0 will  increase  - 
and  in  the  limit,  when  stability  factor  becomes  infinite,  a approaches  1,  the 
motion  becomes  the  vacuum  nutation,  and  points  L and  merge,  0 ceases  to 

have  meaning.  Analogously,  in  precession  (Fig.  l£b)  it  is  the  quasi-nutation 
which  shrinks  more  rapidly,  the  ratio  0/6  decreases  y and  so  does  the  angle  0 
(and  in  the  limit  it  is  f which  ceases  to  have  meaning). 


Passing  to  the  quantitative  evaluation  of  the  rates  of  shrinkage  on  the 
assumptions  of  linearity,  we  must  now  drop  our  pretence  of  avoiding  calculus: 
we  are  now  concerned  with  the  rates  (i.e,,  derivatives)  more  Intimately,  and 
we  are  parallelling  a more  advanced  mathematical  process  (the  USB  perturbation 
for  the  non-constancy  of  the  coefficients  of  the  differential  equations)  than 
before.  Still,  we  shall  need  not  the  theory  of  the  differential  equations, 
but  only  the  most  elementary  operations  of  the  differential  calculus.  If  we 
let  & be  the  quasi-precession  (6  - 0) , we  can  re-write  (3B)  as 

A/6  - (6  - 0)/6  - 1 - 0/6  - (1  l a)/Z  (3C) 

From  Fig.  1 $ we  observe  that  the  rate  of  change  of  A is  A*  **' 
and  the  rate  of  change  of  6 is  6 - - ytj  » it  .will  be  an  easy  matter  to  relate 
the  change  of  velocity  to  the  rate  of  change;  & of  a.  To  relate  these  rates 
(A  8 and  a)  we  need  only  to  differentiate  (3C): 

A /6  - A&/&2  - * */2  <lU) 

which > we  might  note,  we  do. not  consider  as  a differential  equation.  Substi- 
tuting the  expressions  for  A and  6,  and  further  substituting  0 * Y0/A  end 
(Ae»j/B)0  as  before,  we  have  the  equation  for  yt  from  Which 

Y ■ i o/Zifa/A**!  " « (}V) 

Substituting  Into  (lU1)  the  values  of  0/6  and  A /6  from  (IB)  and  (3C),  and 
simplifying  as  before,  we  have 

T « - <V2*(JWj/2B)(l  * a)  (lhA) 

Obviously,  the  proportional  rates  of  change  (6/6)  are 

- a/ 2a  (UiB) 

for  both  precession  and  nutation;  and  little  algebra  need  be  done  to  show 
that  the  radial  rates  with  respect  to  travel  are 

+ (JD+  Je)(l  -l/a2)/2  (ihC) 

2 

(Where  J0  * gdsinS/u  ),  in  agreement  with  the  linear  theory, 

.**  - 

Of  course,  the  mechanism  of  the  resultant  logarithmic  spirals  is  again 
different  from  that  of  the  spirals  of  (8C)  and  IOC),  It  is  indeed  a tribute 
to  the  elegance  of  the  mathematical  linear  theory  that  three  such  distinct 
physical  processes  result  in  the  same  simple  (exponential)  type  of  nol$tdon. 
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It  may  be  remarked  that  the  subscript  D in  (li±C)  refers  to  the  full  drag 
(rather  than  the  axial  drag  DA),  i.e.,  includes  the  yaw-drag;  while  in  the 
standard  texts  on  linear  theory  is  really  meant  as  i.  more  thorough 

(up  to  a point)  study  of  the  theory  indicates , indeed,  that  it  is  more  nearly 
correct  to  use  (in  our  sense)  than  Jda  - as  the  student  will,  no  doubt, 

have  occasion  to  find  out. 

AXIAL  TORQUE 

It  now  remains  only  to  inspect  the  effect  of  the  torque  (ft)  acting  along 
the  axis  of  the  projectile.  In  the  linear  theory  it  is  the  spin-decelerating 
3 2 

torque,  » I^pd  u v,  and  its  effect  is  rather  completely  overshadowed  by  the 

increase  of  stability  due  to  the  loss  of  velocity.  This  torque  can  be  readily 
generalized,  by  the  Maple -Synge  theory,  to  include  terms  analogous  to  yaw-drag, 
etc.* 

Here,  at  last,  is  one  torque  that  truly  maps  on  our  sketches  as  a point  - 
viz.,  the  point  A^.  It  was  mentioned,  however,  that  in  our  method  torques 

must  be  decomposed  along  and  across  L.  In  the  foregoing  discussion  of  the 
cross-torques  - generally,  foUowing**the  assumptions  of  linearity  - we  had 
neglected  this  subtle  distinction,  treating  torques  normal  to  as  though 

they  were  normal  to  Lt  i.e.,  we  had  in  effect  assumed  cos#  - 1,  or  our  error 

was  of  the  second  order  in  #•  The  component  ofi/j  along  L is  approximately 
same  as  , but  its  component  across  L is  approximately  #,  viz.,  of  the 
first  order  in  0 (Fig.  16).  The  firs?  component,  insofar  as  our  sketches  are 
concerned,  acts  analogously  to  a negative  drag,  i.e.,  by  changing  a (except 
that  the  change  of  o due  to  change  of  must  be  distinguished  from  the  change 

of  a due  to  change  of  u)**;  the  second  component  is  readily  seen  to  be  analogous 
to  the  H torque  (except  for  the  sign  and  a constant  factor) . Combining  the 
effects  of  these  two  components,  it  can  be  readily  shown  that  the  proportional 
radial  rates  with  respect  to  travel  are 

(JA/2k|)(i  1/a  * l/o2),  (ISO) 

in  agreement  with  the  linear  theory. 

Again,  certain  refinements  by  the  recission  of  the  assumptions  of  linearity 
are  possible,  if  not  profitable.  In  particular,  it  will  be  readily  seen  that 
there  are  some  effects  due  to  the  components  of  cross-torques  along  L. 


* Cf.  also  BEL  668  and  BRL  882. 

A 

Not  to  be  confused  with  the  fact  that  the  arrow  of  <H  on  Fig*  16  is  in  the 
direction  of  positive  axial  drag. 
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CONCLUSION 


We  have  thus  accomplished  the  following: 

. 1.  Established  an  interpretation,  and  a visualization,  of  the  veU- 
known  results  of  the  mathematical  linear  theory  of  yawing  motion*  This  vis- 
ualization is  not  limited  to  the  linear  theory,  and  may  be  used  as  a starting 
point  for  the  more  advanced  problems  of  ballistics, 

. 2.  Reproduced  all  results  of  the  linear  theory  without  resorting 
to  calculus  (saving  that  excellent  tool  for  the  more  worthwhile  problems  of 
ballistics). 

3*  Reviewed  the  theory  of  yawing  motion  in  the  light  of  the  latest 
experimental  information  (cf .,  e.g.,  the  importance  of  the  non-linearities 
of  T in  comparison  with  XF). 

Outlined  (crudely)  effects  of  some  of  the  non-linearities  of 

ballistics. 

5*  Drawn  attention  to  a certain  inherent  weakness  of  the  present- 
day  mathematical  approach,  viz.,  the  sharp  split  between  the  simple  linear, 
and  generally-unsolved  non-linear,  cases.  It  is  hoped  that  some  day  someone 
will  introduce  some  radical  solution  of  this  difficulty;  it  might  take  the 
form,  for  instance,  of  a special-purpose  analog  differential  analyzer  that 
has,  so  to  say,  the  spherical  trigonometry  built  into  it. 

The  justice  of  the  probable  criticisms,  that  (a)  we  could  not  have  done 
this  work  had  the  mathematical  approach  not  blazed  the  trail,  and  (b)  it  would 
be  simpler  to  learn  calculus,  - is  so  obvious  that  it  requires  no  comment. 

Incidentally,  our  approach  does  seem  to  have  biased  the  trail  in  two  more 
advanced  problems*  the  effect  of  asymmetry  of  the  shell,  (as  mentioned  on 
page  13  of  BRIM  685)  and  the  effect  of  the  liquidity  of  the  filler*  However, 
both  of  those  problems  are  outside  the  scope  of  this  paper. 
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REMARKS  OH  NOTATIONS 


Effort  was  made  to  adhere  to  familiar  notation,  even  at  the  cost  of 
sane  minor  inconsistencies. 

Underlined  letters  stand  for  three-dimensional  vectors.  Otherwise 
letters  are  used  for:  scalar  quantities  (e.g.,  magnitude  of  the  vector); 
to  designate  a point  on  the  sketch;  and  on  a few  occasions,  explained  in  the 
context,  for  complex  quantities. 

Subscript  1 is  used  to  designate  a vector  of  unit  magnitude;  an  exception 
is  the  axial  spin  where  this  subscript  is  used  (after  EEL  hh6)  for  the 

axial  component  of  the  angular  velocity  the  cross -component  is  considered 
either  as  the  vector  T),  or  as  the  scdlsjr'7)  , but  is  visualised  as  an  arrow 
on  the  sketch,  viz.,  essentially  as  the  complex  quantity  f of  BRL  UU6. 

While  the  dot  over  a letter  indicates,  as  usual,  the  derivative  of  a 
quantity  with  respect  to  time,  these  derivatives  (spoken  of  as  "rates*)  are 
introduced  each  time  as  separate  symbols . 

Subscripts  n and  p refer  to  nutation  and  precession. 

Consistent  units  are  assumed  throughout.  All  angles  are  in  radians 
(though  often  units  of  angles  do  not  matter). 

m ■ mass  of  shell;  d - caliber;  u ■ velocity;  g ■ gravity;  8 ■ angle 

of  inclination  of  trajectory,  reckoned -positive  upward;  p«  density  of  the  air; 

A and  5 are  axial  and  transverse  moments  of  inertia,  and  k and  k are  the 

a 

corresponding  radii  of  gyration  in  calibers;  v ■ spin  in  radians  per  caliber 
of  travel;  s - stability  factor,  as  defined  by  (3*)* 

o - / 1 - 1/s,  the  very  convenient  and  familiar  (though  still  nameless) 
parameter  of  exterior  ballistics,  particularly  in  the  linear  theory;  its 
physical  significance  is  clear  fi*om  (3B). 

■ direction  of  trajectory  (velocity  of  the  c.g.  of  shell);  A^  - 

direction  of  shell's  axis;  L ■ angular  momentum  of  the  shell;  G » total  torque 
acting  upon  the  shell.  Letter  0 on  some  sketches  denotes  shell's  center  ef 
gravity. 

6 « yaw  (angle  from  ^ to  A^);  0 « quasi-nutation  (angle  from  1 to  A^; 

also  denoted  by  QN)j  A - quasi-precession  (angle  from  to  L;  also  denoted 

by  QP);  Slf  ■ angular  velocity  of  quasi -nutation  with  respect  to  the  instantan- 
eous position  of  L,  of.  (2);  f(6,0)  • the  trigonometric  function  defined  by 
(U)  and  relevant  For  instantaneously  circular  yawing  motion:  cf . Fig  6;  y 
and  p are, auxiliary  angles,  as  denoted  on  the  sketches. 
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M * overturning  moment;  H ■ damping  moment;  T ■ Magnus  moment;  XT  ■ 
cross-Magnus  moment;  N » Normal  force;  F ■ Magnus  force;  5 ■ Synge  force; 

XF  ® cross-Magnus  force;  D - drag;  DA  • axial  drag.  Italics  A and  «£  are 
used  to  distinguish  the  axial  torque  and  lift  from  the  moment  of  inertia 
and  angular  momentum;  but  little  confusion  would  arise  if  A and  L are  used, 
as  is  customary.  After  Synge,  the  words  "moment*  and  "torque*  are  used 
interchangeably,  as  are  *angular  velocity*  and  *spin*.  K's  are  defined  in 

O O 9 

BBIi  Ui6;  similarly,  J ■ KpdJ/m;  J„  ■ gdsin0/u  and  Jg  ■ gdcos6/u  . As  is 
customary,  coefficient  p.  is  associated  with  M;  to  complete  the  symmetry 
(and  in  minor  variance  from  Fowler  and  EEL  66U),  h is  analogously  associated 
with  H,  t with  T,  and  / with  L. 

Points  of  compass  (N,  NE,  E,  SE,  S,  SW,  W,  NW)  are  used  to  indicate 
the  direction  of  the  yawing  motion,  as  viewed  on  Fig.  U,  qualitatively. 
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